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ABSTRACT 

A  discussion  of  the  physical  nature  of  liquid-film  cooling  is 
presented  that  gives  special  emphasis  to  those  phenomena  occurring  at 
the  gas-liquid  interface.  The  interfacial  mass  and  heat  balances  are 
discussed  in  detail,  together  with  the  phenomena  of  film  instability 
and  the  entrainment  of  liquid  by  the  gas  stream. 

A  review  is  presented  of  the  pertinent  heat-mass  transfer 
analyses  for  the  wall  region  wetted  by  the  liquid  film  The 
limitations  of  those  analyses  are  noted,  and  an  alternate  correlation 
procedure  is  suggested  that  is  applicable  to  the  case  where  the  gas 
flow  is  compressible  and  subject  to  the  influence  of  a  streamwise 
pressure  gradient  In  order  to  implement  that  correlation  procedure, 
a  method  is  developed  whereby  the  temperature  at  the  gas-liquid  inter¬ 
face  can  be  approximated  analytically.  Calculated  results  are  presented 
that  are  noted  to  agree  favorably  with  the  limited  experimental  data. 

A  heat  transfer  analysis  is  presented  for  the  wall  region  downstream 
of  the  liquid  film.  That  analysis  is  based  on  a  relatively  simple  flow 
model,  and  is  applicable  to  the  case  where  the  gas  flow  is  compressible 
and  subject  to  the  influence  of  a  streamwise  pressure  gradient.  Two 
wall  conditions  are  considered:  (a)  the  case  where  the  wall  is  adiabatic; 
and  (b)  the  case  where  the  wall  is  cooled  externally. 

Conclusions  and  recommendations  for  further  experimental  research 


are  presented. 
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1.  INTRODUCTION 
1.1  General  Discussion 

Liquid-film  cooling  refers  to  the  injection  of  a  thin  continuous 
liquid  film  onto  a  given  surface  for  the  purpose  of  protecting  that 
wetted  surface  from  thermal  damage  from  a  hot  gas  flowing  past  the 
film.  I 

The  motivation  for  developing  liquid-film-cooling  technology  for 
rocket  motors  has  been  the  development  of  nuclear  rocket  engines,  the 
use  of  higher  energy  chemical  propellants,  and  the  trend  toward  higher 
combustion  chamber  pressures.  These  newer  developments  are  character¬ 
ized  by  a  large  increase  in  the  heat  flux  from  the  working  fluid  to 
the  walls  of  the  combustion  chamber  and  the  exhaust  nozzle,  such  that 
the  regenerative  cooling  of  those  walls  becomes  marginal  (1) (2) (3) . 

A  need,  therefore,  exists  for  developing  a  method  of  cooling  which 
will  supplement  regenerative  cooling.  While  porous-wall  cooling  tech¬ 
niques  ( transpi ration  cooling  which  utilizes  gaseous  coolants  and  sweat 
cooling  which  utilizes  liquid  coolants)  are  more  effective  than  liquid- 
film  cooling  from  the  standpoint  of  the  degree  of  thermal  protection 
realized  per  pound  of  coolant  consumed,  liquid-film  cooling  still 
appears  to  be  quite  attractive  for  the  following  reasons:  (a)  clogging 
of  the  porous  wall  is  still  a  problem  and  leads  to  localized  "burnout," 
(b)  a  reduction  in  the  structural  integrity  of  large  components  results 
^Numbers  in  parentheses  refer  to  references  listed  in  the  Bibliography. 
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when  they  are  constructed  from  porous  materials,  and  (c)  the  fabrica¬ 
tion  problem  is  simpler  in  the  case  of  a  liquid-film-cooled  system. 

Considerable  attention  has  also  been  devoted  to  studies  con¬ 
cerned  with  employing  a  gaseous  film  coolant  Gaseous-film  cooling 
has  the  same  advantage  of  simplicity  as  does  liquid-film-cooling,  but 
the  absence  of  a  phase  change  during  the  cooling  process  severely 
limits  the  cooling  capacity  of  gaseous-film  cooling  in  comparison  to 
that  of  liquid-film  cooling 

1  2  Objectives  of  the  Present  Investigation 
The  objectives  of  the  subject  investigation  were  threefold: 

(a)  to  describe  the  physical  nature  of  liquid-film  cooling; 

(b)  to  review  the  theoretical  heat-mass  transfer  correlations  for  the 

1 1 qut d- f 1 1m-cooled  region,  and  to  develop  an  alternate  correlation  for 
that  region;  and 

(c)  to  develop  a  heat  transfer  correlation  for  the  region  downstream 
of  the  terminus  of  the  liquid  film 

The  correlation  procedures  developed  herein  are  generalized  so 
as  to  be  applicable  to  the  flow  of  a  compressible  gas  in  the  presence 
of  a  streamwise  pressure  gradient  Moreover,  the  developments  consider 
boundary  layer  flow  rather  than  fully  developed  pipe  flow,  as  have 
previous  analyses  The  developments  are  limited,  however,  to  the  case 
where  the  gas  flow  is  fully  turbulent  and  the  liquid-film  coolant  is 
non-reactive 

Assessment  of  the  validity  of  the  correlation  procedures  developed 
herem  awaits  the  existence  of  accurate  experimental  data.  Such  experi 
mental  work  is  being  conducted  currently  at  the  Jet  Propulsion  Center, 
Purdue  University,  Lafayette,  Indiana. 
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2.  THE  PHYSICAL  NATURE  OF  LIQUID-FILM  COOLING 

2  1  Descri p tio n  _o ■ f_ Liqu id- f  i  1  m  Cooli ng 
Figure  2-1  presents  a  schematic  diagram  illustrating  the  per¬ 
tinent  aspects  of  the  liquid-film  cooling  process.  There  are  two 
regions  of  the  wall  that  are  of  principal  interest:  (a)  the  liquid- 
film-cooled  region,  and  (b)  the  gas-vapor-coded  region  The 
phenomena  which  characterize  each  of  those  regions  are  discussed 
in  the  following  sections 

2.1  1  The  Liquid-Film-Cooled  Region 
Figure  2-2  illustrates  a  physical  boundary  layer  model  for  the 
liquid-film-cooled  region  The  phenomena  which  relate  to  the  develop¬ 
ment  of  that  boundary  layer  region  are  discussed  in  the  following 
paragraphs 

Ihe  Nature  of  Film  Development  At  station  x  =  xQ,  in  Fig  2-1, 
the  liquid  coolant  is  injected  onto  the  solid  surface  and  establishes 
a  spammse  cortmuous  liquid  film  The  l.quid  film  interacts  with 
and  is  accelerated  by  the  hot  qas  stream  flowing  past  it  The  velocity 
distribution  throughout  the  liquid  film  at  any  streamwise  location  is 
determined  by  the  solution  nf  tne  momentum  equation  for  the  two- 
dimensional  flew  of  an  incompressible  rluia;  the  latter  can  be  written 
in  the  followina  form  (4)  Thus 


’■fc* 


mi*. — - 
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The  acceleration  terms  in  Eq.  2-1  are  of  considerable  importance  at 
and  near  the  plane  where  the  liquid  is  injected.  However,  for  the 

largest  portion  of  the  liquid  flow  regime,  because  of  the  extreme 

2 

thinness  of  the  liquid  film,  the  acceleration  terms  and  the  pressure 
gradient  term  in  Eq.  2-1  can  be  considered  of  negligible  importance 
in  comparison  to  the  viscous  term.  Hence,  if  ts  and  tw  denote  the 
shear  stress  at  the  gas-liquid  interface  and  at  the  liquid-solid 
interface,  respectively,  then  one  may  write 


t  h  P]  (v  +  vt)i  constant  •  t$  -  tw 


(2-2) 


where  v  is  the  molecular  contribution  and  vt  is  the  turbulent  or  eddy 
contribution  to  the  "effective"  kinematic  viscosity.  From  Eq.  2-2  it 
follows  that  the  velocity  distribution  for  the  liquid  film  is  given 
by 


u 


y 


r 

i 

o 


dy 


(2-3) 


For  the  special  case  where  the  film  is  entirely  laminar  (vt  =  0), 
Eq.  2-3  reduces  to 


(2-4) 


Refer  to  Appendix  A  for  a  definition  of  the  notation. 

_3 

Typical  film  thicknesses  are  of  the  order  of  1  to  5  X  10  in. 
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The  problem  of  expressing  the  parameter  (v  ♦  vt)j  as  a  function 
of  the  normal  coordinate  y,  so  that  the  integral  term  in  Eq  2-3 
can  be  evaluated  for  the  case  of  a  turbulent  film,  can  be  circumvented 
by  assuming  that  the  universal  "law  of  the  wall"  is  valid  for  the  lquid- 
flow  regime.  cigure  2-3  presents  the  "law  of  the  wall"  in  graphical 
form;  the  dimensionless  parameters  u+  and  y+  presented  in  the  figure 
are  defined  as  follows: 

u+ 

(  VD;ii 

+  y  (  t  /o )  2 
y  = - 

V 

Figure  2-3  represents,  therefore,  a  means  by  which  the  velocity  profile 
in  the  liquid  film  can  be  prescribed  once  the  value  of  the  interfacial 
shear  stress  ’  is  determined;  the  latter  must  result  from  the  solution 
of  the  momentum  equation  for  the  gas-flow  regime. 

The  mass  flow  rate  of  liquid  per  unit  spanwise  length,  denoted 
by  m' ,  and  the  thickness  of  the  liquid  film,  denoted  by  6^,  are  related 
due  to  the  continuity  condition  Hence, 

6f 

m'  =  p|  ;  u  dy  (2-7) 

o 

For  the  case  where  the  film  is  totally  laminar,  Eqs  2-4  and  2-7 
can  be  combined  to  yield  the  following  expression  for  the  thickness 


(2-5) 


(2-6) 
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FIG.  2-3.  THE  UNIVERSAL  "LAW  OF  THE  WALL"  FOR  THE  TURBULENT 
BOUNDARY  LAYER  ON  A  SMOOTH  SURFACE 
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of  the  liquid  film: 


2m 1  ** 

Vvi 


(2-8) 


To  consider  the  relationship  between  m'  and  6^  for  the  case  where  the 
liquid  film  in  turbulent,  it  is  convenient  to  rewrite  Eq.  2-7  in  the 
following  non-dimensional  form.  Thus 


t  +  .  + 

/  u  dy 


(2-9) 


where  the  parameter  W  is  termed  the  dimensionless  liquid  flow  rate, 
and  is  defined  by 


,+  r  m* 


(2-10) 


and  where  6^  is  the  dimensionless  film  thickness  defined  by 


6*  (  T./Pl) 


(2-11) 


The  parameters  W4  and  6^  can  be  related  by  employing  the  "law  of  the 
wall,"  presented  in  Fig.  2-3,  and  graphically  integrating  Eq.2-9. 


- .  siaSms'y .  a. w 
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Figure  2-4  presents  W+  as  a  function  of  6^,  as  determined  by 
Kinney,  et  al.  (5).  The  figure  illustrates,  as  does  Eq.  2-8,  that 
the  thickness  of  the  liquid  film  Is  proportional  to  the  liquid  flow 
rate  m*  and  the  kinematic  viscosity  and  is  Inversely  proportional 
to  the  interfacial  shear  stress  t$.  The  significance  of  the  para¬ 
meter  W+  in  determining  the  physical  characteristics  of  the  flow  of 
a  thin  liquid  film  is  discussed  in  Section  2.4. 

Energy  Transfer  to  the  Liquid  Film.  Energy  is  transferred  to 
the  liquid  film  from  the  "hot"  gas  stream.  A  small  percentage  of 
that  energy  is  utilized  for  increasing  the  sensible  enthalpy  of  the 
liquid  in  the  film.  A  significant  percentage  of  that  energy, 
however,  causes  the  evaporation  of  liquid  which  was  raised  to  its 
saturation  state  while  flowing  downstream. 

If  the  solid  surface  wetted  by  the  film  is  thermally  Insulated 
from  the  external  environment  (i.e.,  adiabatic),  then,  after  that 
surface  is  heated  to  an  equilibrium  temperature  which  is  approximately 
equal  to  that  of  the  liquid  film,  all  of  the  incident  energy  not 
utilized  for  increasing  the  sensible  enthalpy  of  the  liquid  is  utilized 
to  evaporate  it  On  the  other  hand.  If  the  solid  surface  is  externally 
cooled,  then  a  portion  of  the  energy  transferred  to  the  liquid  film 
will  be  transferred  through  the  liquid  film  and  through  the  solid 
surface  to  the  external  coolant.  Thus,  with  external-regenerative 
cooling  of  the  solid  surface,  longer  film-cooled  lengths  will  be 
realized  for  a  given  rate  of  coolant  Injection  than  is  the  case  for 
an  adiabatic  solid  surface. 
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Heat  Transfer  "Blocking.11  The  vapor  produced  by  vaporizing  the 
liquid  film  coolant  is  convected  and  diffused  Into  the  boundary  layer 
of  the  hot  gas  wetting  the  liquid  film  and  establishes  in  turn  a  con¬ 
centration  (or  partial -pressure)  boundary  layer.  In  the  process  of 
being  transported  away  from  the  gas-liquid  interface,  the  vapor 
serves  an  important  function,  referred  to  as  the  "blocking"  of  the 
incident  heat  flux.  In  other  words,  as  the  vapor  moves  away 
from  the  gas-liquid  interface  it  "absorbs"  sensible  heat,  so  that 
it  thus  represent^  a  counter-convective  heat  transfer  mechanism. 

Fi 1m  termination.  At  station  x  *  Xp  in  Fig.  2-1,  the  liquid 
film  terminates,  due  to  evaporation  and  possible  entrainment  of  liquid 
droplets  into  the  gas  stream.  An  undesirably  high  degree  of  droplet 
entrainment  is  known  to  arise  when  large  wave  crests,  termed  distur¬ 
bances,  occur  on  the  surface  of  the  liquid  film.  The  phenomenon 
of  large  disturbances  on  the  surface  of  the  film  is  termed  f i lm 
instability  and  is  discussed  in  Section  2.4. 

Coupling  of  the  Transport  Processes  and  the  Boundary  Conditions. 
The  interdependence  of  the  heat  and  mass  transfer  processes  and  the 
related  boundary  conditions  for  the  evaporation  process  is  illustrated 
by  the  following  points: 

(a)  The  rate  at  which  energy  is  transferred  to  the  liquid  film  and 
the  physical  properties  of  the  liquid  (particularly  Its  latent  heat 
of  vaporization)  determine  the  rate  at  which  vapor  Is  "Injected"  into 
the  gaseous  boundary  layer. 
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(b)  The  rate  of  "vapor  injection"  determines  the  degree  of  the 
aforementioned  "blockage"  phenomenon  and  thus  the  net  rate  of 
energy  transfer  to  the  liquid  film. 

(c)  The  concentration  of  the  vapor  at  the  interface  is  dependent 
on  the  rate  at  which  vapor  is  introduced  into  the  boundary-layer 
region  (see  Section  2.2.1,  Eq.  2-23). 

(d)  The  temperature  at  the  interface  is  determined  by  the  con¬ 
centration  of  the  vapor  at  the  interface  (or,  perhaps  more  cor- 
recfly,  the  partial  pressure  of  the  vapor  at  the  interface). 

(e)  The  rate  at  which  energy  is  transferred  to  the  liquid  film  is 
dependent  on  the  thermal  potential  which  exists  between  the  main 
stream  and  the  gas-liquid  interface. 

The  physics  of  the  evaporation  process  thus  require  that  for  a 
given  set  of  bulk  flow  conditions  (i.e.,  working  fluids,  main-stream 
temperature, pressure,  Reynolds  number,  Mach  number,  etc.),  unique 
values  for  each  of  the  following  must  result:  the  rate  at  which  energy 
is  transferred  to  the  liquid  film;  the  rate  at  which  liquid  is 
evaporated;  and  the  concentration  of  the  vapor  and  the  temperature 
at  the  gas-liquid  interface. 

2.1.2  The  Gas -Vapor-Cooled  Region* 

Figure  2-5  illustrates  schematically  what  is  defined  as  the  gas- 

vapor-cooled  region,  together  with  representative  profiles  of  the 

wall  temperature  for  the  cases  where  the  wall  Is  adiabatic  and  is 

1  The  terminology  "vapor-cooled  region"  Introduced  in  Refs  (6)  and  (7) 
would  appear  to  be  a  misnomer  In  that  the  quantity  ot  relatively 
cool  gas  which  exists  throughout  the  boundary  layer  at  x*x,  contributes 
substantially  to  the  gas-vapor-cooled  length.  1 
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GAS  STREAM  — 


LIQUID  FILM 


GAS-VAPOR-COOLED  LENGTH 
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ADIABAT'C  OR  E XTERNALLY- COOLED  WALL 


t. 


FIG.  2-5.  A  SCHEMATIC  REPRESENTATION  OF  THE 
PERTINENT  WALL  TEMPERATURE  PROFILES 
FOR  THE  GAS-VAPOR-COOLED  REGION 
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externally  cooled  The  lenqth  of  the  gar  -vapor-cooled  region,  given 
by  x~-x,  in  the  figure,  'S  defined  as  that  in  which  the  temperature 

C  i 

of  an  adiabati;  surface,  denotea  by  t  ,  increases  from  the  tempera¬ 
ture  of  the  liquid  film  at  x  -  x. ,  denoted  by  t$  to  the  recovery 
temperature  for  the  main  stream 

Ii  the  surface  downstream  from  the  end  cf  the  liquid  film  is 
cooled  externally,  then  a  wall  temperature  distribution  such  as  that 
denoted  by  t  in  Fig  2-5  results.  In  general,  at  any  station 
x  >Xj,  the  rate  of  heat  transfer  to  the  external  coolant  is  given  by 


where  h'  is  the  appropriate  heat  transfer  film  coefficient  to  be 
w 

employed  with  the  thermal  potential  (tflw  -  t  )  in  the  calculation  of 
the  wall  heat  flux  qw  To  either  design  an  effective  regenerate vely- 
cooled  downstream  region  or  determine  the  optimum  location  for  intro¬ 
ducing  a  second  liquid-film-cooled  region,  it  must  be  possible  to 

determine  the  streamwise  distribution  of  q  for  a  prescribed  distribo- 

w 

tion  of  t  Hence,  the  basic  problem  is  that  of  determining  the 
streamwise  distribution  of  the  adiabatic-wall  temperature,  together 
with  the  heat  transfer  film  coefficient  An  analysis  for  the  gas- 
vapor-cooled  region  is  presented  in  Action  4  which  considers  those 
two  points. 
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2,2  Conservation  of  Mass  and  Energy  at  the- Gas-Liquid  Interface 
The  conservation  of  mass  and  energy  across  the  gas-liquid  inter¬ 
face  is  considered  in  Sections  2=2.1  and  2.2.2,  respectively.  Con¬ 
sideration  is  given  only  to  the  case  where  the  liquid  film  is  non¬ 
reactive  and  is  stable,  free  from  the  consequences  of  large  disturbances 
(see  Section  2.4).  It  is  not  assumed  that  the  film  is  either  smooth 
(free  from  waviness)  or  laminar.  Moreover,  mass  diffusion  in  the 
boundary  layer  is  treated  as  a  binary  process;  i.e.,  the  evaporated 
vapor  and  the  gas  entrained  from  the  main  stream  into  the  boundary 
layer  are  considered  as  the  two  counter-diffusing  species. 

2,2.1  Interfacial  Mass  Balance 

Figure  2-6  illustrates  a  general  interfacial  element  for  mass- 
transfer,  together  with  the  several  mass  fluxes  which  are  discussed 
in  the  following  paragraphs.  The  S-  and  L-surfaces  presented  in 
Fig.  2-6  are  defined  as  being  placed  in  the  gas-vapor  phase  and  in 
the  liquid  phase,  respectively,  at  an  infinitesimal  distance  from 
the  interface.  Additional  definitions  and  subscript  notation  needed 
to  clarify  Fig  2-6  are  as  follows: 

Definitions: 

(pv>,  *  Vv  *  °v 

lov)9  *  Vs  =  '^d.S1 


(2-12) 

(2-13) 


(PV)  =  (Pg  +  Py)  V 


r  i 

=  the  partial  density  for  the  "i“  species;  i=g,v 

i 


(2-14) 
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(pv)  (p  v  ) 

v  S  v  d,v  S 


He  H  c  ^Vl  ^V.  )e 

S  v,S  g  S  g  d,g  S 


(yov) 


L 


FIG.  2-6.  A  MASS-TRANSFER  INTERFACE  ELEMENT 
FOR  A  NON-REACTIVE,  STABLE  LIQUID 
FILM 
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-  the  partial  pressure  for  the  "i"  species;  i=g,v 

o=  oy  +  r  the  density  for  the  gas-vapor  mixture 
v^  -  v  ♦  vi  d  ;  the  "total"  velocity  normal  to  the  S-  and  L- 
suf faces  for  the  "i“  spec'es;  i-g,v 
v  -  the  convective  velocity  normal  to  the  S-  and  L-surfaces 
v-  d  =  the  diffusional  velocity  normal  to  the  S-  and  L-surfaces 
for  the  "i "  species  :  i rg ,v 

Subscripts 

v  -  vapor 

g  -  gas  entrained  from  the  main  stream 
L  -  evaluated  at  the  L-surface 
S  -  evaluated  at  the  S-surface 

d  -  diffusional  component  of  the  subscripted  parameter 

A  tenn  having  neither  "v,"  "g"  nor  "L"  as  a  subscript  refers  to  the 
gas-vapor  mi  xtu-'-e 

F’qure  2-6  ’Hustrates  that  when  the  liquid  coolant  is  evaporating 
into  a  medium  other  than  its  own  vapor,  mass  transport  is  due  to  both 
diffusion  and  convection.  The  relationships  governing  the  inter¬ 
dependence  or  the  different  convective  and  diffusive  mass  fluxes 
presented  in  Fig  2-b  are  derived  in  the  following  paragraphs 
Picks  law  relates  the  mass  transfer  by  diffusion  m  the  y 
direction  to  the  concentration  gradient  in  the  y  direction,  and  is 
expressed  by 

dC. 

(ovd)i  -d  (D*Dt)  '  {i  v,g)  (2-15) 
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where 

o 

C  -  -  -  the  concentration  of  the  "i"  species  (i  =  g,v). 


(D  D t )  1  the  "effective"  diffusion  coefficient 


In  the  case  of  a  binary  medium,  to  which  the  subject  discussion  is 
restricted,  the  concentration  gradients  aC^/ and  3C^/ ay  are 
related  by 


dC 

v 


(2-16) 


Equations  2-15  and  2-16  show  that  at  any  station  y  the  rates  of  mass 
diffusion  of  the  vapor,  'py^,  and  the  gas,  (ov^)  ,  are  related  by 


Since  the  interface  is  impermeable  to  the  gas,  the  diffusive  flux 
of  gas  towards  the  liquid  surface  is  exactly  equal  to  the  convective 
flux  of  qas  away  from  the  liquid  surface;  i  e  , the  net  transport  of 
gas  at  the  surface  of  the  liquid  film  is  zero.*  thus 


( c  v  . )  c 

d  g,S 


I  . 

u 


v). 
<j  5 


(2-lb) 


Hence 


d  ,g,S 


=  -  vr 


(2-19) 


Note,  however,  that  the  concentration  of  the  gas  is  always  greater 
than  zero  at  the  interface 


?n 

Equations  2-17  and  2-19  can  be  combined  to  yield 

K’v.S  ■  O  -  Cv,S»»v>S  (2-20) 

Equations  2-12  and  2-20  can  be  combined  to  yield 

(pv)s  =  (Pv)v>s  (2-21) 

Hence,  Eqs.  2-20  and  2-21  show  that  at  the  interface,  the  diffusive 
rate  of  vapor  transport,  Uv,j)v  $♦  and  the  total  rate  of  vapor 
transport,  (ov)  are  related  by 

V  ,j 

<«d>»,S  *  <l  -  Cv.S»»»)v,S  <2-221 

It  follows  from  the  conservation  of  the  evaporating  species 
across  the  interface  and  from  Equations  2-21,  2-22,  and  2-15,  that 

p(D+DJ  aC  ~ 

<»*>,. s‘  <«»»L  *  <«>s  ■  •  -fir-sr  {z-n) 

-  v  -  s 

Equations  2-15,  2-17,  2-19,  2-22,  and  2-23  describe,  for  the 
process  of  evaporation,  the  pertinent  relationships  between  the  several 
mass  fluxes  presented  in  Fig  2-6  Those  equations  also  can  be 
written  in  terms  of  the  parti  a1  pressure  of  the  vapor  Consider 
Eq  2-15:  from  the  definitions  of  the  concentration  of  the  vapor, 

C  ,  and  the  partial  pressure  of  the  vapor,  Py,  ft  follows  that 

aC  ,  pR  n  ap  C  . D 

_ v  _  3_  _v__\  _  R _  rv  v  aR 

ay  '  ay  p  R/  ‘  p  \  ay  R  ay 


(2-24) 
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2.2,2  Interfacial  Energy  Balance 
Figure  2-7  illustrates  a  general  energy-transfer  interface 
element  placed  at  the  surface  of  a  stable,  non-reactive  liquid  film. 
The  different  energy  fluxes  Involved  are  defined  as  follows: 

(1)  '  (k  |y)^=  the  rate  at  which  energy  is  transferred  across 
the  S-surface  by  conduction 

(2)  qL  =  (k  — )L  =  the  rate  at  which  energy  is  transferred  across 
the  L-surface  by  conduction. 


(3)  (ov)  «.  h  <.  =  the  rate  at  which  energy  is  convected  across 
the  S-surface  by  the  evaporated  vapor. 

(4)  (ov)  hL  =  the  rate  at  which  energy  is  convected  across  the  L- 
surface  by  the  evaporating  liquid. 

(5)  q  <.  *  the  rate  at  which  energy  is  transferred  across  the 

r 


(6) 


S-surface  by  radiation 

q  -  the  rate  at  which  energy  is  transferred  across  the 

r  ,l 


L-surface  by  radiation 

US  VS 

(7)  (KE)S  (°v)s  (-j-  +  -y)  =  the  rate  at  which  kinetic  energy 

associated  with  the  gas-vapor  phase  crosses  the  S-surface. 

UL  VL 

(8)  (KE)l  (pv)l  {-£  ♦  -j-)  *  the  rate  at  which  kinetic  energy 
associated  with  the  liquid  phase  crosses  the  L-surface. 


Since  the  net  transport  of  the  gas  at  the  surface  of  the  liquid 

film  is  zero,*  only  the  transport  of  enthalpy  associated  with  the 

evaporated  vapor  (see  Item  (3)  above)  needs  be  considered  on  the 

gas-vapor  side  of  the  interface  It  is  desirable,  therefore,  that 
*See  discussion  leading  to  Eq  2-18. 


(KE)s  is 


q  ip\i)  c  h 
r,S  r  v, S  v,S 


S -SURFACE 


i  GAS- VAPOR  PHASE 

LIQUID  PHASE 


INTERFACE 


L- SURFACE 


till 

q  q  (KE)  ip\i)  h 
L  r,L  Lr  LL 


FIG.  2-7.  AN  ENERGY-TRANSFER  INTERFACE  ELEMENT 
FOR  A  STABLE,  NON-REACTIVE  LIQUID  FILM 
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the  vapor  have  a  relatively  large  specific  heat  so  that  the  vapor 
will  cause  a  large  amount  of  thermal  "blockage." 

In  view  of  Eq.  2-23,  the  energy  balance  for  the  Interface 
element  illustrated  in  Fig.  2-7  can  be  written  as  follows: 

«s  *  \  -!iKE)S  •  <KE>l|*  (»'>S(hv,S-hL)-(Vs-<>r,L>  <2-30) 

I  II  III  IV  V 

Term  IV  in  Eq.  2-30  can  be  rewritten  in  the  form 

(pv)s  &Hy 

where  AHy  is  the  latent  heat  of  vaporization  for  the  liquid  coolant 
which  corresponds  to  the  saturation  conditions  prescribed  by  the 
partial  pressure  of  the  vapor  at  the  interface.  If  the  partial 
pressure  of  the  vapor  at  the  interface  is  equal  to  or  greater  than 
the  critical  pressure  for  the  liquid  coolant,  then  term  IV  in 
Eq  2-30  is  identically  zero.  In  most  instances,  term  III  In 
Eq.  2-30  can  be  considered  negligible  in  comparison  to  a  term 
such  as  term  I  Thus,  neglecting  term  III,  Eq.  2-30  can  be  rewritten 
to  yield 

qL  *  V,L  '  «S  *  VS  '  ‘*»>S  lHv  <2-31) 

Equation  2-31  is  the  general  energy  balance  for  conditions  at  the 
interface.  To  obtain  a  better  understanding  of  what  Is  Involved 
in  evaluating  the  left-hand  side  of  Eq.  2-31,  the  energy  balance  for 
the  entire  liquid  film  must  be  considered,  and  that  energy  balance 
is  developed  in  the  following  section. 
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2.3.1  General  Energy  Balance 
Figure  2-8  illustrates  diagrammatically  a  control  volume 
placed  within  a  stable,  non-reactive  liquid  film.  The  following 
energy  fluxes  may  be  defined:  ? 

(1)  m;  h1  o  =  the  rate  of  enthalpy  transfer  into  the  control 
volume  per  unit  spanwise  length  due  to  the  continuous  in¬ 
jection  of  film  coolant. 

(2)  qw  =  (k  —  )w  =  the  rate  of  conduction  of  heat  from  the 
liquid  film  to  the  wall. 

(3)  q  =  the  rate  of  radiant  energy  transfer  to  the  wall. 

r  ,w 

The  rate  of  injection  of  liquid  into  the  control  volume 
illustrated  in  Fig.  2-8  per  unit  spanwise  length  is  denoted  by  m^. 
From  the  conservation  of  liquid  film  coolant  it  follows  that 


/  (ov).  dx 

xo 


(2-32) 


Employing  Eq.  2-32,  the  energy  balance  for  the  control  volume 
presented  in  Fig.  2-8  is  given  by 


'  (q,  *  q  . )  dx 

Xq  xhL  1  r ,  L ' 


kq  L,‘,v)S(hL'hl,o)  +  qw  +  ^r,w  <2-33) 


Equation  2-33  is  the  general  energy  balance  for  the  liquid  film. 
For  the  special  case  where  the  summation  of  the  terms  under 
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NON-REACTIVE  LIQUID  FILM 
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each  of  the  integral  signs  in  Eq.  2-33  is  independent  of  x,  the 
integral  signs  can  be  removed  and  one  obtains  a  local  energy  balance 
for  the  liquid  film.  Furthermore,  the  relationship  obtained  after 
removing  the  integral  signs  can  be  combined  with  Eq.  2-31,  obtaining 
thereby  a  loca7  energy  balance  at  the  interface.  Thus 

%  ♦  Vs  =  "w  +  Vw  +  (<*v>S(hL-  hl  ,o)H°v)S  aHv  <2-34> 

In  many  liquid-film  cooling  applications  of  interest,  however, 
the  terms  under  the  integral  signs  in  Eq.  2-33  are  not  independent 
of  the  streamwise  coordinate  x.  Of  particular  interest  is  the 
problem  of  liquid-film  cooling  an  exhaust  nozzle  wherein  severe 
pressure  gradients  are  common,  and  wherein  the  various  terms  in 
Eq.  2-33  are  most  certainly  functions  of  x.  In  such  a  situation 
it  is  difficult  to  appraise  the  error  involved  when  it  is  assumed 
that  Eq.  2-34  is  the  local  energy  balance  at  the  interface.  In 
particular,  it  is  difficult  to  appraise  the  error  due  to  assuming 
that  the  portion  of  the  incident  heat  flux  utilized  for  increasing 
the  sensible  enthalpy  of  the  liquid  is  equal  to  (dv)^  (h^  -  h^  Q) , 
as  given  by  Eq.  2-34.  Since  Eq.  2-34  is  much  more  convenient  to 
employ,  however,  than  is  Eq.  2-31,  further  consideration  of  the 
problem  is  worthwhile.  A  method  for  experimentally  investigating 
that  problem  is  suggested  in  Section  3.3.1. 
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2  3,2  Simplified  Energy  Balances  Employed  by  Previous  Investigators 

Knuth  (8)  simplified  the  general  interfacial  energy  balance  (see 
Eq  2-31)  by  assuming,  in  essence,  that 

qr,S  ‘  qL  r  qr,L 

In  other  words,  Knuth  equated  the  radiant  energy  incident  on  the  sur¬ 
face  of  the  liquid  film  to  the  energy  utilized  for  increasing  the 
sensible  enthalpy  of  the  liquid  coolant  plus  the  energy  transferred 
to  the  solid  surface  wetted  by  the  film.  (In  Knuth 's  experimental 
investigation,  the  wetted  surface  was  approximately  adiabatic  ) 

Knuth,  therefore,  employed  the  following  interfacial  energy  balance 
Thus 

qs  »  (ov)$  AHv  (2-35) 

Crocco  (9)  utilized  the  same  interfacial  energy  balance  as 
Knuth,  but  Crocco 's  reasoning  was  somewhat  different.  He  assumed 
that  the  rate  of  radiant  energy  trarsfer  to  the  liquid  film  was 
negligible  and  that  the  wetted  surf  ce  was  adiabatic  Crocco 
further  assumed  that  in  a  so-called  "first  period  of  existence," 
all  of  the  film  coolant  was  heated  to  equilibrium  saturation 
conditions,  and  that  for  the  major  portion  of  the  liquid- flow 
regime,  all  of  the  energy  incident  on  the  liquid  film  was  utilized 
for  evaporating  it.  Moreover,  the  analysis  was  concerned  only  with 
that  major  portion  of  the  liquid-flow  regime 
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Graham  (10),  Sellers  (11),  Kinney  (5),  and  Emmons  (6)  employed 
an  interfacial  energy  balance  which  is  a  simplification  of  Eq.  2-34. 
First,  it  was  assumed  that  the  radiant  energy  transfer  was  negligible. 
Furthermore,  in  the  experimental  investigations  conducted  by  Graham, 
Kinney,  and  Sellers,  and  in  a  portion  of  the  experimental  Investigation 
conducted  by  Enmons,  the  wetted  surface  was  essentially  adiabatic. 

Thus,  the  interfacial  energy  balance  employed  by  those  investigators  for 
the  case  where  the  wetted  surface  was  adiabatic  is 

qs  =  (pv)<.  (hL*hl,o^  +  (pv)S  A  Hv  (2-36) 

2.4  Gas-Liquid  Interactions 

The  specification  of  the  physical  characteristics  of  the  gas- 
liquid  interface  is  of  primary  importance  in  the  analysis  of  liquid- 
film  cooling.  Considerable  experimental  data  are  available  that 
indicate  that  the  interfacial  structure  has  a  considerable  influence 
on  the  flow  rate  of  film  coolant  required  for  establishing  a 
specified  liquid-film-cooled  length  (5) (8) ( 10) ( 11 ) .  Of  particular 
significance  is  the  phenomenon  of  film  Instability  which  causes  the 
entrainment  of  liquid  coolant  droplets  by  the  flowing  gas  stream. 

The  physical  nature  of  the  gas-liquid  Interface  is  discussed  in 
Section  2.4.1  and  the  phenomena  of  film  instability  and  entrain¬ 
ment  are  discussed  in  Section  2.4.2. 

2.4.1  Interfacial  Structures 

There  appear  to  be  three  distinct  interfacial  structures 
associated  with  the  process  of  liquid-film  cooling  that  are  of 
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interest.  Figure  2-9  presents  photographs  of  each  of  those  inter¬ 
facial  structures.  The  photographs  were  taken  during  certain  pre¬ 
liminary  investigations  conducted  at  the  Jet  Propulsion  Center, 

Purdue  University  (12).  The  following  features  of  the  basic  flow 
system  and  of  the  photographic  apparatus  are  noteworthy:  (a)  the 
gaseous  and  liquid  phases  were  air  and  deionized  water,  respectively; 
(b)  the  gaseous  flow  was  turbulent;  (c)  a  35  mm  camera  equipped  with 
a  close-copy  lens  was  employed  in  conjunction  with  a  stroboscopic 
light  source,  which  had  a  flash  time  of  2  microseconds,  so  that  the 
gas  side  of  the  two-phase  interface  could  be  photographed;  and  (d) 
the  optical  line  of  sight  was  at  an  angle  of  10°  with  respect  to  the 
wetted  test  surface  which  was  a  highly  polished,  flat,  stainless 
steel  plate  having  a  spanwise  dimension  of  4  ins.. 

The  photographs  presented  in  Fig.  2-9  were  taken  when  the  static 
pressure  and  the  main-stream  temperature  had  values  of  220  psia  and 
1000  R,  respectively.  Flow  parameters  which  were  not  common  to  all 
of  the  photographs  are  specified  in  the  appropriate  subtitles  of 
Fig.  2-9;  they  are  the  main- stream  velocity,  ue>  the  position  down¬ 
stream  of  the  point  of  liquid-film  coolant  injection  at  which  the 
photographs  were  taken,  x  -  xQ,  and  the  rate  of  liquid-film  coolant 
injection  per  inch  of  injector  slot  width,  m^. 

The  interfacial  structure  presented  in  Fig.  2-9a  is  characterized 
by  its  extreme  smoothness.  In  that  particular  photograph  there  is 
shown  a  H  in.  diameter  cylindrical  plug  which  protrudes  approximately 
0.001  in.  above  the  surface  of  the  polished  stainless  steel  surface 
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o  i,n 

X-SCOle:  I _ I _ I 


a.  smooth  (ue=88fps  ,  x-x0=9ins  ,  m^sO-OOl  Ib/in-soc) 


b."pebbled"  l  ue=50fps  ,  x-x0=5ins  ,  mJfO-002  Ib/in-sec) 


c.Hp«k)bled"  (uesl30fps,  x-x^Sins  ,  rr!^  O-OOI  Ib/in-sec) 


d.  unstable  (  ue=  130  fps  ,x-x0=  5in$  ,  0  0045  Ib/in-sec) 

FIG.  2-9.  CHARACTERISTIC  INTERFACIAL  STRUCTURES 
FOR  LIQUID-FILM  COOLING 
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The  fortuitous  presence  of  that  surface  irregularity  serves  to 
indicate  that  a  liquid  film  which  has  a  smooth  interfacial 
structure  is,  in  general,  very  thin  and  sensitive  to  relatively 
small  surface  roughnesses.  The  smooth  character  of  the  inter¬ 
facial  structure  is  apparently  due  to  the  fact  that  the  liquid  film 
is  so  thin  and  viscous  that  any  disturbances  imparted  to  it  as  a 
consequence  of  turbulent  velocity  fluctuations  in  the  adjacent  gas 
phase  are  completely  suppressed  or  "damped  out."  Because  of  the 
aforementioned  sensitivity  of  a  thin  film  to  surface  roughnesses, 
it  is  highly  improbable  that  the  liquid  films  employed  for  liquid- 
film  cooling  the  walls  of  a  rocket  thrust  chamber  will  have  an 
extremely  smooth  interfacial  structure. 

The  interfacial  structure  presented  in  Figs.  2-9b  and  2-9c 
is  quite  common  and  is  characterized  by  its  "pebbled"  appearance 
The  small  disturbances  which  are  evident  on  the  surface  of  the  liquid 
film  are  apparently  a  consequence  of  turbulent  velocity  fluctuations 
in  the  gas  phase  flowing  past  the  liquid  film  (13).  The  liquid  film 
is  sufficiently  thin  and  viscous  that  small  disturbances  which  appear 
on  its  surface  do  not  amplify  and  cause  film  instability.  (It  should 
not  be  implied,  however,  that  there  may  not  exist  turbulent  stresses 
in  the  liquid  film.)  The  difference  in  the  scale  of  the  interfacial 
disturbances  for  the  two  "pebbled"  interfaces  presented  in  Fig  2-9 
is  due  primarily  to  a  corresponding  difference  in  the  scale  of  the 
turbulence  in  the  adjacent  flowing  gas;  the  turbulence  scale  of  the 
gas  phase  decreases  with  increasing  gas  mass  flow  rate.  That  relation¬ 
ship  between  the  gas  mass  flow  rate  and  the  nature  of  the  "pebbled" 
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interfacial  structure  has  been  observed  previously  (5) (8) (13X 14) . 

The  interfacial  structure  presented  in  Fig.  2-9d  defines  the 
phenomenon  of  fi lm  instabi 1 i ty;  the  latter  is  characterized  by 
the  existence  of  large-scale  disturbances  which  are  turbulent  in 
nature,  and  which  are  superimposed  on  a  "pebbled”  interfacial 
structure  The  phenomenon  of  film  instability  and  the  related 
phenomenon  of  liquid  entrainment  are  discussed  in  the  following 
section. 

2.4.2  Film  Instability  and  Entrainment 

A  large  degree  of  liquid  entrainment  is  known  to  occur  when 
large-scale  disturbances  appear  on  the  surface  of  a  liquid  film 
(15).  Such  entrainment  reduces  the  effectiveness  of  the  liquid- 
film  cooling  process  in  either  or  both  of  two  ways: 

1.  If  the  entrained  coolant  is  convected  out  of  the  region  which 
is  to  be  liquid-film  cooled  before  it  is  evaporated,  then  a  cor¬ 
responding  loss  in  liquid-film  cooling  effectiveness  arises  due  to 
the  incomplete  utilization  of  the  available  heat  of  vaporization 
associated  with  the  flow  of  the  liquid-film  coolant. 

2.  If  the  entrained  liquid  droplets  are  evaporated  immediately  upon 
leaving  the  interface,  then  the  corresponding  multi-fold  increase 

in  the  specific  volume  of  the  coolant  upon  evaporation  will  have  a 
disturbing  effect  on  the  gaseous  boundary  layer  structure  which, 
in  turn,  might  tend  to  promote  the  degree  of  turbulent  heat,  mass, 
and  momentum  transport  near  the  gas-liquid  interface.  Luikov  (16, 
17)  has  termed  such  a  phenomenon  volumetric  boiling. 
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Thus,  a  problem  of  primary  importance  is  that  concerned  with 
the  specification  of  (a)  the  physical  conditions  under  which  small- 
scale  interfacial  disturbances  will  amplify  and  cause  film 
instability,  and  (b)  the  degree  of  and  the  effects  of  liquid 
coolant  entrainment  once  transition  has  occurred.  Intuitively, 
parameters  which  may  be  of  importance  in  a  formulation  of  the 
solution  to  the  above  problem  are  as  follows: 

(i)  the  iriterfacial  shear  stress,  rg,  which  represents  the  degree 
to  which  the  gaseous  flow  regime  tends  to  deform  the  liquid  film; 

(ii)  the  dynamic  viscosity  of  the  liquid,  p-j ,  which  is  a  measure  of 
the  resistance  of  the  liquid  film  to  deformation; 

(iii)  the  thickness  of  the  liquid  film,  6^,  which  represents  the 
extent  to  which  the  gas-liquid  interface  experiences  the  stabilizing 
influence  of  the  wetted  solid  surface; 

( i v)  the  kinematic  viscosity  of  the  liquid,  ,  which  is  a  measure 
of  the  damping  characteristics  of  the  liquid  film; 

(v)  the  interfacial  density  ratio,  ps/pj»  which  represents  the 
relative  inertia  characteristics  of  the  two  phases,  and  is  a  measure 
of  the  degree  to  which  a  turbulent  velocity  fluctuation  in  the  gaseous 
phase  will  influence  the  liquid-flow  regime; 

(vi)  the  dynamic  viscosity  ratio  p<j/  which  is  a  measure  of  the 
discontinuity  in  the  velocity  profile  at  the  gas-liquid  interface; 
and 

(vii)  the  surface  tension  of  the  liquid  phase. 

References  (5) (13) ( 14) ( 15)  have  experimentally  investigated  the 
phenomenon  of  film  instability.  Subsequent  attempts  to  correlate  the 
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inception  of  film  instability  have  employed  the  dimensionless  liquid- 
flow  rate  W+  (see  Eq.  2-10).  It  is  of  interest  to  note  that  W+  com¬ 
bines  the  individual  influences  of  the  parameters  (1 )-(1v)  listed 
above  into  a  single  parameter. 

The  motivation  for  attempting  to  correlate  the  inception  of  film 
instability  in  terms  of  the  single  parameter  W+  arises  from  the 
following  facts  which  characterize  single  phase  turbulent  flow: 

(a)  a  universal  relationship  between  u+  and  y+  (the  "law-of-the-wall") 
characterizes  the  "inner"  region  (y+  <  1000)  of  the  wall  turbulent 
region1;  and  (b)  the  distance  normal  to  the  wall  at  which  turbulent 
forces  are  of  the  same  order  of  importance  as  the  viscous  forces 
is  defined  approximately  by  15<y+  <20  and  90  <  W+ <  150  (5) ( 18) . 

Thus,  the  transition  state  for  the  liquid-flow  regime  would  be 
characterized  by  corresponding  limits  of  6^  and  W+  if  the  following 
basic  hypothesis  is  physically  correct:  (1)  parameters  (v)-(vii), 
presented  in  the  foregoing,  do  not  appreciably  affect  the  transi¬ 
tion  characteristics  of  the  liquid-flow  regime;  (2)  the  liquid-flow 
regime  obeys  the  universal  "law-of-the-wall"  as  established  for 
single  phase  flow;  and  (3)  transition  occurs  when  the  turbulent 
forces  in  the  region  of  the  liquid  film  close  to  the  gas-liquid 
interface  are  of  the  same  order  of  importance  as  the  viscous  forces. 

While  a  large  percentage  of  the  experimental  data  tends  to 
substantiate  the  foregoing  simple  hypothesis  (5) ( 15) ,  a  limited 
amount  of  data  has  been  obtained  which  indicates  that  the  simple 
hypothesis  may  not  give  a  complete  characterization  of  the  phenomena. 
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For  example,  experimental  observations  indicate  that  the  critical 
value  of  W+  increases  slightly  with  decreasing  surface  tension  (5), 
and  that  the  parameter  u^/u^  may  also  be  significant  in  defining 
the  transition  state  ( 13) ( 14) .  Furthermore,  at  relatively  low  values 
of  the  gas  stream  Reynolds  number *(<10^)  it  is  known  that  the 
critical  value  of  W+  can  increase  significantly  (19);  i.e.,  the 
transition  state  is  a  function  of  both  the  Reynolds  number  for  the 
gas  stream  and  W+.  It  would  appear,  however,  that  at  values  of  the 
gas  stream  Reynolds  number  which  are  most  commonly  associated  with 
the  gas  flow  in  a  rocket  thrust  chamber  ( >10^ ) ,  the  dependence  of 
the  critical  value  of  W+  on  the  Reynolds  number  for  the  gas  stream 
is  negligible  (5 ) ( 14) ( 19) . 

For  the  practical  design  of  a  liquid-film-cooled  system  which 

5 

is  characterized  by  gas-stream  Reynolds  numbers  greater  than  10  , 
the  following  generalized  results  given  by  Kinney,  et  al,  (5)  appear 
to  be  adequate  until  a  more  complete  theory  is  formulated  and  sub¬ 
stantially  verified: 

1.  In  the  situation  where  the  evaporation  rate  is  negligible,  and 
wherein  m'  is  thus  independent  of  the  streamwise  coordinate  x,  the 
inception  of  film  instability  is  approximately  defined  by  W*  =  90 
where 


The  viscosity  in  Eq.  2-37  is  evaluated  at  a  streamwise-mean 
liquid-film  temperature. 

1  The  Reynolds  number  referred  to  on  this  page  is  R.,  defined  on 
page  57.  a  . 
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2.  In  the  case  of  liquid-film  cooling,  wherein  the  value  of  W+ 
varies  from  a  nonzero  value  at  x  =  xQ  to  zero  at  x  =x1  ,  the 
onset  of  a  high  degree  of  liquid  entrainment--which  accompanies  a 
corresponding  decrease  in  liquid-film  cooling  effectiveness--!' s 
approximately  defined  by  W*  =  360. 

Graham  (10,20)  experimentally  investigated  the  relationship 
between  the  effectiveness  of  liquid-film  cooling  and  the  phenomenon 
of  liquid  entrainment.  He  introduced  a  stability  effectiveness  e$ 
defined  by 


e 


s 


m‘ 


o^i 

o,a 


(2-38) 


where 

m'o  j  =  the  ideal  rate  of  coolant  consumption  (zero  entrainment)  . 

m'o  a  ■  the  actual  rate  of  coolant  consumption. 

Graham  correlated  es  in  terms  of  W*  utilizing  the  experimental  data 
of  References  (5)  and  (10);  the  resultant  correlation  for  smooth- 
tube  flow  is  presented  in  Fig.  2-10.  It  should  be  noted  that  the 
correlation  was  developed  utilizing  the  assumption  that  the  heat 
transfer  film  coefficient  h^  is  independent  of  W*.  There  appear 
to  exist,  however,  considerable  experimental  data  which  indicate 
that  h^  is  significantly  dependent  on  W*;  the  latter  phenomenon  is 


discussed  in  Section  3.2. 
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FIG.  2-10.  THE  CORRELATION  DUE  TO  GRAHAMOO) 
FOR  THE  STABILITY  EFFECTIVENESS  €. 
AS  A  FUNCTION  OF  THE  LIQUID  FLOW 
RATE  PARAMETER  W+ 
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The  stability  effectiveness  ■  represents  a  measure  of  an 
effect  of  the  entrainment  phenomenon  rather  than  a  direct  measure 
of  the  actual  quantity  of  coolant  which  is  entrained  by  a  flowing 
gas  stream.  The  latter  problem  has  been  experimental ly  investigated 
by  Hewitt,  et  al .  (21)  for  the  case  of  annular  two  phase  flow,  in 
that  subject  investigation  the  quantily  of  entrained  liquid  was 
determined  by  employing  a  sampling  probe  which  traversed  the  central 
core  of  the  two-phase  flow  field.  Ii  was  observed  that  below  a 
critical  value  of  liquid  flow  rate,  the  entrainment  was  negligible; 
the  critical  state  was  approximately  defined  by  W*  =  90  (15). 

As  the  liquid  flow  rate  increased  above  the  critical  value,  the 
trainment  of  liquid  into  the  central  gas  core  increased 
correspondingly.  It  was  observed,  however,  that  the  degree  of 
entrainment  increased  with  increasing  rate  of  gas  flow  as  well  as 
with  increasing  rate  of  liquid  flow.  Furthermore,  additional 
experiments  conducted  at  different  system  pressures  indicated  that 
the  degree  of  entrainment  increased  with  decreased  pressure,  the  mass 
flow  rate  of  gas  being  held  constant.  Those  two  facts  would  indicate 
that  the  interfacial  shear  stress  r  ,  in  addition  to  the  dimensionless 
liquid  flow  rate  W*,  might  be  of  significance  in  determining  the  actual 
quantity  of  coolant  which  is  entrained  by  a  flowing  gas  stream. 

The  values  cf  the  experimental  parameters  investigated  by  Hewitt, 
et  al.  (21)  were  significantly  different  from  those  which  would 
normally  characterize  a.i  application  of  liquid-film  cooling.  In 
particular,  the  temperature  and  pressure  employed  were  essentially 
ambient  and  the  evaporative  rate  of  mass  transfer  was  zero  It  is 
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impossible,  therefore,  to  determine  from  those  data  the  relation¬ 
ship  which  exists  between  the  stability  effectiveness  cg 
and  the  actual  quantity  of  coolant  which  is  entrained  by  a  flowing 
gas  stream.  It  would  appear  that  a  need  exists  for  additional 
experimental  investigation  of  the  phenomenon  of  liquid  entrainment 
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3.  ANALYZING  THE  L I QUI D- F I LM- COOLED  REGION 

A  review  of  the  previous  analyses  for  liquid-film  cooling  is 
presented  in  Section  3.1,  and  a  review  of  some  pertinent  experimental 
data  is  presented  in  Section  3  2.  The  latter  serves  to  illustrate  the 
limitations  of  those  analyses  that  are  reviewed,  and  provides  the  moti¬ 
vation  for  the  alternate  correlation  procedure  that  is  presented  in 
Section  3  3  for  heat-mass  transfer  data  that  are  obtained  in  liquid- 
film  cooling  experiments 

3 . 1  A  Review  of  the  Previous  Analyses 
Table  1  compares  the  analyses  due  to  Knuth  (8),  Seller.;  (11) (22), 
Graham  (10),  Ranm'e  (23) \  and  Emmons  (6).  Column  I  of  that  table  shows 
that  similar  sets  of  equations  are  employed  in  the  analyses  due  to 
Knuth,  Sellers,  Graham,  and  Ranm'e;  consequently,  those  analyses  are 
considered  concurrently  in  Section  3.1  1.  Because  the  analysis  due 
to  Emmon^  i.,  trie  most  recent  analysis  developed,  and  also  since  it 
differs  r«.idamental  ly  from  the  earlier  analyses,  that  analysis  is 
considered  separately  in  Section  312 

The  fundamental  steps  that  are  common  to  each  of  the  subject 
analyses  are  outlined  briefly  below: 

Step  1.  Elementary  boundary-1 ayer  theory  is  employed  to  analyze  the 
gaseous-flow  regime  above  the  liquid  film  In  the  analyses  due  to  Knuth, 
*  As  modified  for  liquid- film  cooling  by  Sloop  (24). 
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Rannie,  Sellers,  and  Graham,  the  Prar.dtl-Taylor  theory  (18)  is  extended 

to  include  the  effects  of  mass  transfer,  in  the  analysis  due  to  Emmons, 

the  theory  developed  by  Turcotte  (25)  for  the  case  of  the  transpired 

turbulent  boundary  layer  on  a  flat  plate  is  extended  to  the  case  of 

liquid-film  cool’ng  The  result  of  chose  developments  is  basically  an 

i  2 

analytical  expression  for  the  following  non-dimensional  parameter  ’  . 
Thus 


(ov)s  (te-ts) 


(3-1) 


where 


3 

q<-  =  the  rate  at  which  energy  is  transferred  across  the  S-surfare 

by  conduction. 

( o v ) c  =  the  rate  at  which  mass  is  transferred  across  the  S-surface 
by  convection 

=  the  temperature  at  the  edge  of  the  boundary  layer,  or  on  the 
axis  of  the  duct. 

t  =  the  tempei iture  at  the  gas-liquid  interface, 

c  =  the  specific  heat  for  the  vapor. 
pv 


Step  2.  Another  expression  is  developed  for  the  parameter  6  by  consider- 

4 

i.ig  the  conservation  of  energy  across  the  gas-liquid  interface 

*  The  parameter  .<  was  denoted  by  H  in  the  original  analyses  due  to  Sellers, 
and  is  r:luted  to  the  Stanton  number  for  heat  transfer,  defined  by 
Eq  3-25 

'  See  Table  1  rrj  uir.n  1 1 1 

^  See  Section  2  2  1  for  the  definition  of  the  S-surface. 

4 

See  Table  1,  column  V 
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That  expression  for  g  combined  with  the  expression  for  B  obtained  in 
Step  1  yields,  in  principle,  a  single  equation  in  the  single  unknown 
(pv)s. 

Step  3.  Finally,  from  the  conservation  of  mass  for  the  liquid-film 
coolant,  the  coolant  flow  rate  m^  is  related  to  the  liquid-film-cooled 
length  L  by 

m'  =  (pv)$  L  ( 3-2 ) 1 

Equation  3-2,  which  relates  m^  and  L,  is  what  is  desired  as  the  result 
of  each  analysis. 

In  completing  Step  1,  which  constitutes  the  major  part  of  each 
analysis,  there  arises  more  unknown  quantities  than  there  are  inde¬ 
pendent  equations.  To  circumvent  that  difficulty,  one  of  the  following 

2 

methods  is  employed: 

(a)  Empirical  formulae,  based  on  zero  mass  transfer,  are  employed  to 
obtain  approximate  values  for  the  unknown  quantities  needed  for  completing 
the  solution  The  appropriateness  of  the  assumptions  employed  is  deter¬ 
mined  by  comparing  the  analytical  predictions  with  experimental  data  (8) 
(23). 

(b)  Some  of  the  unknown  quantities  needed  for  completing  the  solution 
are  approximated  by  means  of  empirical  formulae,  and  the  remainder  are 
determined  experimentally.  An  obvious  disadvantage  of  the  latter  method 

1  Equation  3-2  follows  from  Eq.  2-32  and  Assumption  (l.f)  of  the  analyti¬ 
cal  model  that  is  presented  in  the  following  section. 

2 

See  Table  1,  column  IV. 
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is  that  the  final  analytical  expression  cannot  be  applied  with  confi¬ 
dence  to  situations  that  differ  significantly  from  that  for  which  the 
unknowns  were  determined.  Sellers  determined  two  unknowns  and  Graham 
and  Emmons  each  determined  one  unknown  from  experimental  data. 

3.1.1  The  Analyses  Due  to  Rannie,  Sellers,  Graham,  and  Knuth 
Discussed  in  the  present  section  are  the  fundamental  developments 
that  relate  to  the  completion  of  Step  1  for  the  subject  analyses. 

The  Physical  and  Analytical  Models.  Figure  3-1  illustrates 
schematically  a  physical  model  characteristic  of  those  employed  by 
References  (8)(11)( 22) ( 10) ( 23) .  The  assumptions  common  to  the  different 
analytical  models  are  as  follows: 

(1) 1  The  partial  derivatives  of  the  following  parameters  with  respect  to 
the  streamwise  coordinate  x  are  zero: 

(l.a)  the  streamwise  component  of  velocity,  u; 

(l.b)  the  static  temperature  of  the  gas  stream,  t; 

(l.c)  the  partial  pressure  of  the  vapor,  py; 

(Id)  the  interfacial  shear  stress,  ts; 

(l.e)  the  rate  at  which  heat  is  transferred  across  the  S-surface 
by  conduction,  and 

(l.f)  the  rate  at  which  mass  is  transferred  across  the  S-surface 
by  convection  (pv),.. 

(2)  The  gaseous  flow  regime  is  subdivided  into  an  "effective"  laminar 
sublayer,  bounded  by  y  =  0  and  y  *  6,  and  a  turbulent  core,  bounded  by 
y  =  6  and  y  *  A. 

*  Assumption  (1)  is  equivalent  to  either  the  assumption  of  fully- 
developed-pipe  flow  or  the  assumption  of  Couette  flow. 
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FIG.  3-1.  A  SCHEMATIC  REPRESENTATION  OF  THE 
PHYSICAL  MODEL  EMPLOYED  IN  THE 
SUBJECT  ANALYSES 
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(3)  The  wetted  surface  is  adiabatic. 

(4)  The  liquid-film  coolant  is  non-reactive,  and  the  interfacial  sur¬ 
face  is  stable. 

(5)  The  flowing  gas  is  incompressible,  and  its  properties  may  be 
evaluated  at  a  suitable  reference  temperature. 

(6)  The  effects  of  mass  transfer  upon  heat  transfer  and  momentum 
transfer  are  important  only  in  the  laminar  sublayer  (except  in  Knuth's 
analysis  where  said  effects  are  also  considered  in  the  turbulent  core). 

(7)  The  turbulent  Prandtl  number  and  Schmidt  number,  defined  by 


and  Sct 


respectively,  are  unity. 


(8)  The  surface  velocity  of  the  liquid  film  Is  negligible  compared  to 
the  main-stream  velocity  up. 

(9)  Across  the  "thin"  laminar  sublayer  it  can  be  assumed  that 


(9- a) 

(pv) 

-DRp  dp 

=  tt—j --  - — -i  -  ~~  -  constant  * 
RV(P-PV)  dy 

(,v)s 

(3-3) 

(9  b) 

du 

U  dy  " 

(pv)su  =  constant  *  x$ 

(3-4) 

(9.  r) 

- 

dy 

(pv),  c  t  *  constant  *  q,  - 
3  Py  5 

(pw,s  V  *s 

( 3-5 ) 1 

It  is  argued  in  Ref.  (26),  pp.  10-12,  that  the  specific  heat  for 
the  vapor  in  Eq.  3-5  should  be  replaced  by  the  specific  heat  for 
the  gas-vapor  mixture.  However,  that  argument  Is  incorrect  due 
to  the  fact  that  the  diffusive  mode  of  mass  transport  was  not 
considered  (see  the  developments  leading  to  Eqs.  2-18  and  2-30). 
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The  Resistance  of  the  Laminar  Sublayer  to  Momentum  Transfer.  I t 
follows  ‘hat  by  integrating  Eq.  3-4  and  substituting  the  limits  u  =  0 
at  y  =  0  and  u  =  ufi  at  y  =  6,  one  obtains  an  equation  for  the  resistance 
of  the  laminar  sublayer  to  momentum  transfer-  In  dimensionless  form, 
that  equation  is 


o  u 

e 


(pu)e 

2  (pvjs 


(pv)s6 
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The  Resistance  of  the  Laminar  Sublayer  to  Heat  Transfer.  It 
follows  that  by  integrating  Eq.  3-5  and  introducing  the  limits  t  =  t, 
at  y  =  0  and  t  =  tfi  at  y  =  6,  one  obtains  an  equation  for  the  resis¬ 
tance  of  the  laminar  sublayer  to  heat  transfer.  In  non-dimensional 
form  that  equation  is 


t,  -  t  , 

t - t-  =  \  (exp  (Pr 

*e  ls  8 


Ms 

v 


6)  -1) 


(3-7) 


The  Resistance  of  the  Laminar  Sublayer  to  Mass  Transfer.  Of  the 
subject  analyses,  only  that  analysis  due  to  Knuth  includes  a  speci¬ 
fication  of  the  vapor-partial-pressure  field  (or  its  equivalent, 
the  vapor-concentration  field).  In  general,  the  specification  of  that 
field  may  be  desirable  for  one  or  more  of  the  following  reasons:  (a) 
the  evaluation  of  the  physical  and  transport  properties;  (b)  the  evalua¬ 
tion  of  the  energy  released  by  chemical  reactions  in  the  case  where  the 
liquid-film  coolant  is  reactive;  and  (c)  the  determination  of  the 


temperature  at  the  gas-liquid  interface--which  is  uniquely  related  to 
the  partial  pressure  of  the  vapor  at  the  interface.  If  that  tempera¬ 
ture  must  be  known  with  a  high  degree  of  accuracy,  the  approximation* 
introduced  by  Graham,  Sellers,  and  Emmons  is  inappropriate.  Knuth 
considered  the  problem  of  specifying  the  partial-pressure  field  for 
the  latter  reason. 

For  completeness,  the  specification  of  the  partial  pressure  dis¬ 
tribution  is  considered  herein.  Integrating  Eq.  3-3  and  imposing  the 
limits  p  =  p  at  y  =  0  and  p  =  p  at  y  =  6,  the  equation  is 

V  V  ,S  V  V  $ 6 

obtained  for  the  resistance  of  the  laminar  sublayer  to  mass  transport, 
which  in  non-dimensional  form  is  given  by 
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i he  Relative  Resistance  of  the  Turbulent  Core  to  Mass,  Momentum , 
and  Heat  Transfer  (the  Reynolds  Analogy).  A  logical  extension  of  the 
Reynolds  analogy  (18) (27)  to  heat,  mass,  and  momentum  transfer  in  the 
cur'uu ierit  CGiC  .jaicuub-i  iuv,!  regiuu.  i:> 


(mass  flux)  _  (momentum  flux)  _  (heat  flux)  _ 

(mass  flux)  _  (momentum  flux)  _  (heat  flux)  _  g 

y  Jr  y 

The  above  relationship  can  be  written  more  explicitly  in  two  ways: 

(a)  based  on  the  diffusive  rates  of  transport  only,  and  (b)  based  on 

the  total  rates  of  transport  (diffusion  plus  convection).  Rannie, 

Sellers,  and  Graham  employed  method  (a);  Knuth  employed  method  (b). 

*  The  approximation  being  that  the  temperature  of  the  liquid  film 
equals  the  saturation  temperature  uf  the  film  coolant  which  cor¬ 
responds  to  the  static  pressure  p. 


The  Reynolds  analogy  which  considers  only  the  diffusive  mode 


of  transport  can  be  written  as  follows,  Thus 


-D, 


Rip 

v  r 


Tpv); 


du 

ut  dy  = 
*6 


(3-9) 


where 

=  the  rate  at  which  energy  is  transferred  across  the  surface 
y  =  6  by  conduction 

(pv)f  =  the  rate  at  which  mass  is  transferred  across  the  surface 
y  =  6  by  convection, 

r ,  =  the  shear  stress  at  y  =  6 

6 


Integrating  Eq  3-9,  substituting  the  appropriate  limits  at  y  =  6 
and  at  y  =  a,  and  introducing  Assumption  (7),  yields 
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(b) 


(c) 


The  Reynolds  analogy  that  cinsiders  both  convective  and  diffusive 
rates  of  transport  can  be  written  in  the  following  form  (8): 
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Integrating  Eq.  3-11,  substituting  the  appropriate  limits  at  y  *  6 
and  at  y  =  a,  and  employing  Assumption  (7),  one  obtains  (8) 


P  -  P 


*  (pv), 

-  In  U  +  — £  (u  -  u.)  } 


(3-12) 


-P-  In  {1  + 

?v 


(te  -  V  ‘„v> 


Equations  3-10  and  3-12  are  the  two  ways  in  which  the  Reynolds 
analogy  can  be  employed  to  describe  the  relative  resistance  of  the 
turbulent  core  to  the  transfer  of  heat,  mass,  and  momentum.  It  should 
be  noted  that  Assumption  (9)  (see  page  47),  which  was  employed  to 
derive  the  equations  for  the  laminar  sublayer,  is  a  special  case  of 
the  Reynolds  analogy  that  considers  both  diffusive  and  convective 
transport  (Eq.  3-11)  To  consider  two  of  the  analyses  due  to  Sellers, 
the  Reynolds  analogy  that  considers  only  diffusive  heat  and  momentum 
transport  must  also  be  written  for  the  laminar  sublayer;  i.e.. 


k  dy  _  dt  _  a  constant  across  the. 
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Integrating  Eq.  3-13  across  the  sublayer  and  substituting  the  appro¬ 
priate  limits,  one  obtains* 

1  Sellers  states  that  Eq  3-14  implies  that  the  velocity  and  temperature 
profiles  in  the  laminar  sublayer  are  assumed  to  be  linear  (11,  p.  20); 
however,  since  Eq.  3-13  needs  not  be  so  restricted,  it  would  appear 
that  qualifying  states  t  by  Sellers  is  unnecessary. 
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Comparing  the  Analyses,  Equations  3-3  through  3-8  plus  Eqs.  3-10, 
3-12,  and  3-14  represent  the  relationships  from  which  the  analyses  due 
to  Rannie,  Sellers,  Graham,  and  Knuth  have  been  developed  Table  1, 
p.  42,  indicates  the  specific  equations  that  are  employed  in  each  of 
the  subject  analyses  (column  I),  together  with  the  unknowns  that  are 
eliminated  between  those  equations  (column  II)  and  the  resultant 
analytical  expression  for  the  parameter  8  (column  III)  Also  presented 
in  that  table  are  the  empirical  formulae  and/or  experimental  data  that 
are  employed  in  each  of  the  subject  analyses  (column  IV),  and  the 
interfacial  heat  balance  that  is  employed  (column  V) 


312  The  Analysis  Due  to  Emmons^ 

Physical  and  Analytical  Models.  The  physical  and  analytical 
models  presented  in  Section  3  1  I  are  applicable  to  the  analysis  due 
to  Emmons,  with  the  exception  of  Assumptions  (2)  and  (9),  Assumption  (2) 
should  be  rewritten  as  follows:  The  gaseous-fluw  regime  is  subdivided 
into  a  turbulent  core  and  a  viscous  sublayer  region  wherein  the  following 
eddy-vi scosi ty  relationship  proposed  by  Turcotte  (25)  applies: 
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*  To  improve  the  clarity,  completeness,  and  correctness  of  the  original 
analysis  due  to  Emmons,  the  developments  and  the  final  solution 
presented  herein  differ  from  that  originally  presented  by  Emmons; 
however,  the  analytical  and  physical  models  are  the  same. 
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The  parameter  b  in  Eq.  3-15  must  be  determined  experimentally.  ^ 

Assumption  (9)  should  be  modified  slightly  in  accordance  with  the 

new  definition  of  the  sublayer  region;  i,e.,  the  molecular  thermal 

conductivity  k  and  dynamic  viscosity  u  should  be  replaced  by  their 

"effective"  values,  denoted  by  k  +  kt  and  w  +  ut,  respectively. 

An  additional  significant  assumption  employed  by  Emmons  is  that 

v  c 

the  molecular  Prandtl  number,  defined  by  Pr  =  .  is  equal  to  unity. 


Tg  Related  to  q<.  (The  Reynolds  Analogy).  To  relate  the  Inter¬ 
facial  shear  stress  to  the  heat  flux  q<-,  Emmons  assumed  that  the 
basic  Reynolds  analogy,  written  as  follows,  is  valid  for  the  entire 
gaseous-flow  regime: 
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Integrating  Eq.  3-16  and  substituting  the  appropriate  limits  at  y  *  0 
and  at  y  -  A,  one  obtains 


(3-17) 


Evaluating _ s  To  evaluate  the  Interfacial  shear  stress 

Emmons  employs  a  relationship  that  was  derived  by  Turcotte  (25,  Eq.  15) 
for  the  case  of  the  transpired  turbulent  boundary  layer  on  a  flat  plate. 


*  The  parameter  b  was  introduced  by  Turcotte  Into  the  relationship 
originally  proposed  by  Rannie  (28)  for  the  case  of  zero  mass 
transfer;  the  parameter  b  Accounts  for  the  fact  that  it  is  not 
completely  clear  as  to  whether+T  is  the  correct  shear  stress  to 
employ  in  the  calculation  of  y  . 
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where 

x  =  the  shear  stress  at  the  wall  in  the  absence  of  mass  transfe 
s  ,o 

As  indicated  in  Table  1,  Eqs  3-17  and  3-18  constitute  an  analytical 

.  .  2 

solution  for  the  parameter  8.  Since  Emmons  misquotes  and  misinterprets 
Fq.  3-18,  it  is  of  value  to  outline  Turcotte's  derivation  of  that 
equation.  By  combining  Eqs.  3-4  and  3-15,  Turcotte  obtains  the  follow¬ 
ing  equation  for  the  shear  stress  distribution  in  the  sublayer  region. 
Thus 

T  s  u  (1  +  sinh2  (fj^gg))  *  Ts+  (pv)s  u  (3-19) 


Integrating  Eq.  3-19  and  substituting  the  limit  u  =  0  at  y  =  0,  one 
obtains  the  equation  due  to  Turcotte  for  the  velocity  profile  in  the 
viscous  sublayer.  Thus 
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Combining  Eqs  3-19  and  3-20,  one  obtains  the  following  equation  for 

the  shear  stress  profile  in  the  viscous  sublayer: 

1  Emmons  writes  Eq  3rl8  incorrectly  He  replaces  the  quantity 
( t  /p )*»  by  (t  / cP  (see  Ref  (6),  Eq  B-ll) 
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1  See  Ref  (6),  pp  84,  85,  91 
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iVrotte  then  int'oduces  the  assumpt  on  that  the  turbulent  core  ex- 
terns’  to  the  viscous  sublayer  region  is  unaffected  by  the  mass 
transfer  at  the  wan  Since,  *n  the  absence  of  mass  transfer, 
the  wa’l  region  defined  by  y+  1000  is  a  constant-stress  layer 
( 18) ,  it  fol lows  that 


and  that 
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limit  tanh  ^3 -gg— )  1  1  (3-23) 
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Equations  3-22  and  3-23  thus  reduce  Eq  3-21  to  Eq  3-18. 

Comparison  of  Emmons 1  Analysis  with  the  Analyses  Due  to  Sellers, 
Graham,  Rannie,_and  Knuth  Table  1,  p  42,  compares  the  analysis  due 
to  Emmons  with  those  analyses  d’scussed  in  Section  3  1.1  Column  III 
01  that  table  presents  two  analytical  express’''. is  for  the  paramber  e 
under  the  analysis  due  to  Emmons:  (a)  the  expression  due  to  Emmons 
(6,  Eq  B- 19 ) ;  and  (b)  the  corrected  expression  derived  by  the  writer 

32  A  Review  of  the  Pertinent  Experimental  Data 
The  analyses  reviewed  in  Section  3  l  are  oased  on  the  premise  that 
when  the  liquid  film  is  stable,  the  analytical  end  physical  mode's  fur 
liquid-film  cooling  and  t  ranspi  ration  coolmg  may  br  considered 
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substantially  identical.  That  is,  the  surface  of  the  liquid  film  may 
be  treated  as  a  zero- velocity ,  aerodynamically  smooth  porous  surface 
through  which  the  vapor  is  "injected"  into  the  adjacent  boundary  layer. 
Moreover,  the  concept  of  the  viscous  sublayer  is  introduced,  that 
concept  having  been  employed  with  success  in  the  analysis  of  trans¬ 
piration  cooling.  A  review  of  some  pertinent  experimental  data  is 
presented  here  to  demonstrate  that  there  is  a  serious  question  regarding 
the  validity  of  that  premise.  Discussed  in  Section  3.2.1  are  heat- 
mass  transfer  data  that  indicate  both  the  limitations  of  those  theories 
reviewed  in  Section  3  1  and  the  complex  nature  of  the  phenomena  that 
characterize  the  evaporation  process.  Data  that  present  an  explana¬ 
tion  for  anomalies  that  are  noted  in  the  heat-mass  transfer  data  are 
discussed  in  Section  3.2.2. 

3.2  1  Heat-Mass  Transfer  Data 

The  following  items  are  discussed  In  the  present  section: 

(a)  the  agreement  between  the  heat-mass  transfer  rates  predicted  by 
the  theory  due  to  Knuth  (8)  and  the  data  reported  by  Refs.  (5)  and  (8); 

(b)  the  physical  significance  of  the  correlation  that  was  obtained  by 
Emmons  (6)  for  his  experimental  data;  and 

(c)  the  data  obtained  by  Luikov  (17)  in  an  investigation  of  the  basic 
phenomena  that  characterize  the  evaporation  of  a  liquid  In  the  presence 
of  a  "hot"  gas  stream 

A  Comparison  of  Knuth 's  Theory  with  Experimental  Data.  Knuth 
observed  that  the  average*  heat-mass  transfer  rates  predicted  by  his 
*  Averaged  over  the  liquid- fi Im-cooled  length. 
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theory  were  significantly  smaller  than  those  he  obtained  experimentally, 
and  those  Kinney,  et  al.  (5)  reported;  the  difference  was  attributed  by 
Knuth  to  "entrance  length"  effects.*'  These  were  not  accounted  for  in 
any  of  the  other  analyses  reviewed  in  Section  3.1.  Knuth  employed  an 
empirical  correction  to  bring  the  theoretical  and  experimental  heat- 
mass  transfer  rates  into  agreement;  that  empirical  correction  consisted 
of  multiplying  the  theoretical  heat  and  mass  transfer  rates  by  the 
following  expression.  Thus  \ 

O  ♦irRd*> 

Rd  =  the  Reynolds  number  for  the  gas  stream  based  on  the  diameter  of 
the  tube  and  the  bulk  flow  parameters, 
d  -  the  diameter  of  the  tube. 

L  =  the  liquid-film-cooled  length. 

Typical  values  of  Rd  and  d/L  for  the  data  analyzed  by  Knuth  were  10^ 
and  10_1,  respectively.  Hence,  the  empirical  correction  employed  bv 
Knuth  was  of  the  order  of 

1  +  1/3  10"1  [\06)k  =  2.05 

The  average  heat-mass  transfer  rates  predicted  by  Knuth* s  theory  were, 
therefore,  typically  one  half  of  those  observed  experimentally. 

T~  ‘  "  “ 

*  In  the  experimental  investigations  conducted  by  Knuth  and  Kinney, 
et  al.,  the  velocity  profile  at  the  plane  of  liquid  injection  was 
essentially  fully  developed,  whereas  the  temperature  and  concentra¬ 
tion  profiles  were  flat. 
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and  the  thermal  boundary  conditions  are  such  that  the  velocity  profile 
at  the  entrance  to  the  test  section  is  fully  deveioped  while  the  tempera¬ 
ture  profile  is  flat,  the  empirical  correction  for  the  "entrance  length" 
effects  is  given  by  the  following  expression  (8): 


1+0.1- 


where  L  is  the  length  of  the  test  section  over  which  heat  transfer  occurs. 
Thus,  if  R,  =  10^  and  d/L  =  10"*,  the  ratio  of  the  average  rate  cf  heat 
transfer  over  the  test  section  length  L  to  the  heat  transfer  rate  that 
would  obtain  if  the  temperature  profile  also  was  fully  developed  is 

1  +  10" 1  (106)’4  =1.32 

It  would  appear,  therefore,  that  the  presence  of  the  evaporating  liquid 
film  in  the  experimental  investigations  conducted  by  Knuth  and  Kinney, 
et  al.  resulted  in  a  substantial  increase  in  the  average  heat  transfer 
rate  over  that  which  would  obtain  if  the  mass  transfer  was  zero,  that 
ratio  being  typically  2.05/1.32  =  1.55  Basic  mass  transfer  theory, 
however,  would  predict  that  the  "injection"  of  the  vapor  into  the  gas 
stream  would  act  to  "block"  the  heat  transfer  to  the  liquid  film,  thus 
reducing  the  observed  rates  of  heat  transfer  to  the  film  over  those 
that  would  result  for  the  case  of  zero  mass  transfer;  the  data  for 
transpiration  cooling  presented  in  Fig  3-3  in  Section  332  illustrates 
that  point  It  would  appear,  therefore,  that  a  significant  anamolie 
exists  in  the  data  reported  for  liquid-film  cooling  by  Refs.  (5)  and  (8) 
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The  Physical  Significance  of  the  Correlation  Due  to  Eranons  (6) . 

As  noted  in  Table  1,  the  value  of  the  parameter  b  with  which  Emmons 
was  able  to  correlate  his  experimental  data  was  0.110  However,  on  a 
physical  basis,  Turcotte  (25)  argued  that  b  must  be  bounded  as  follows: 


where,  for  transpiration  cooling, 

( »  /  -  )  ^  >  1 

'  s,o/  s; 

Turcotte  was  able  to  obtain  $atisfuCto«\>  correlation  of  experimental 
data  for  the  transpired  turbulent  boundary  layer  by  employing  the 
median  value  of  b;  thus 


Obviously,  the  value  cf  b  employed  by  Emmons  for  correlating  his 
experimental  data  disagrees  significantly  w’th  the  va^es  of  b  that  are 
physically  justified.  Some  of  the  discrepancy  may  be  due  to  the  in¬ 
correct  quotation  of  Eq  3-18  by  Emmons,  as  discussed  in  Section  3.1  2 
It  would  appear,  however,  that  a  portion  of  the  discrepancy  is  attribut¬ 
able  to  an  incorrect  assumption  in  his  analytical  model  The  discussion 
presented  in  the  following  section,  Section  322,  indicates  that  the 
eddy-viscosi ty  law  given  by  Eq  3-15,  which  forms  the  basis  for  Emmons' 
correlation,  does  not  generally  apply  to  the  gaseous-flow  regime  ad¬ 
jacent  to  the  surface  of  the  liquid  film 
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The  Evaporation  Data  Reporter  r.y  luikov  (17)  The  fcl Icwir.y  poif.U 
concerning  the  experimental  apparatus  and  the  experimental  procedure 
employed  in  the  subject  investigation  are  worthy  of  note: 

(a)  A  dry  body  which  was  internally  cooled  and  a  shallow  pan--having 
similar  overall  dimensions  as  the  dry  body--contai ning  a  liquid  were 
placed  side  by  side  on  the  floor  of  a  large  wind  tunnel. 

(b)  The  liquid  in  the  evaporating  pan  was  conti nuous ly  replenismn. 
so  that  the  free  surface  of  the  liquid  was  always  flush  with  the  top 
of  the  pan. 

(c)  Both  the  free  surface  of  the  liquid  in  the  evaporating  pan  and 
the  surface  of  the  internally  cooled  dry  body  were  fully  instrumented 
with  thermocouples ;  in  all  of  the  experiments,  the  rate  of  coolant  flow 
through  the  dry  body  was  adjusted  so  that  the  temnerafu'-e  of  the  free 
liquid  surface  and  the  surface  of  the  dry  body  were  the  same 

(d)  Four  liquids  were  employed  in  the  investigation:  benzol,  acetone, 
butanol ,  and  water. 

(e)  The  static  pressure  in  the  wind  tunnel  was  essentially  ambient 
and  the  temperature  of  the  main  stream  varied  from  40  C  to  120  C, 
and  the  velocities  from  5  rreters/sec  to  14  meters/sec. 

The  experimental  results  were  as  follows:  in  all  cases,  the 
observed  average  heat  transfer  coefficient  for  the  free  liquid  surface 
was  larger  than  the  average  heat  transfer  coefficient  for  the  dry  body. 
The  ratio  of  the  two  heat  transfer  coefficients  varied  from  1.2  to 
1.6--the  larger  values  corresponding  to  the  larger  ovaporaMo..  rates. 

As  indicated  in  the  foregoing  paragraphs,  that  trend  in  the  experimental 
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data  is  opposed  to  the  trend  which  would  be  predicted  by  mass  transfer 
theory.  Moreover,  the  similarity  in  the  thermal  boundary  condition 
that  was  established  in  the  subject  investigation  for  both  the  dry 
surface  and  the  surface  of  the  evaporating  liquid  discredits  to  a  large 
degree  the  proposal  by  Knuth  that  the  disagreement  between  his  theory 
and  experimental  liquid-film  cooling  date  is  attributable  solely  to 
"entrance  length"  effects. 

32  2  Correlation  of  Velocity  Profile  ta 
In  the  present  section  the  data  reported  by  Refs  (19)  and  (21) 
are  noted;  they  relate  to  the  phenomenon  of  gas-liquid  shear  inter¬ 
actions  in  tne  absence  of  evaporative  mass  transfer  The  experimental 
apparatus  employed  in  each  of  the  subject  investigations  were  signifi¬ 
cantly  different;  for  a  complete  description  of  either  experimental 
apparatus  the  reader  is  referred  to  the  references  The  experimental 
results,  however,  were  in  qualitative  agreement  and  may  be  summarized 
as  follows  When  traverses  were  made  with  a  pitot  probe  across  the  gas- 
flow  regime  adjacent  to  a  "pebbled"  gas-liquid  interface,^  the  resultant 
velocity  profile  data  correlated  in  the  same  manner  as  velocity  profile 
data  for  rough-pipe  flow  It  was  ooserved  that  the  "effective  roughness" 
of  the  liquid  film  surface  increased  with  increasing  values  of  the 
dimensionless  liquid  flow  rate  W+  (defined  by  Eq.  2-10) 

Because  the  viscous  sublayer  affords  the  bulk  of  the  resistance  of 
the  boundary  layer  to  the  transport  of  heat  from  a  hot  gas  stream  to  a 

*  Refe>"  to  Fig  2-9,  presented  in  Section  2-4,  for  a  photographic 
definition  of  the  "pebbled"  interface  structure. 


62 


solid  boundary  (18),  and  since  the  sublayer  region  becomes  less 
influential  and  less  well  defined  as  the  surface  roughness  increases 
(4)  (18),  it  would  appear  that  the  interfacial  phenomenon  discussed 
above  may  account  to  a  large  degree  for  the  anomalies  in  the  heat- 
mass  transfer  data  discussed  in  the  previous  section. 

3.3  An  Alternate  Correlation  for  Heat-Mass  Transfer  Data 
Consider  the  following  recapitulation  of  the  experimental  informa¬ 
tion  presented  in  Sections  3  2  and  2  4  For  the  case  where  the  liquid 
film  is  stable,  the  effectiveness  of  1'qu  d-film  cooling  apparently 
decreases  with  increasing  values  of  the  dimensionless  liquid  flow  rate 
W+,  due  to  a  corresponding  increase  in  the  "effective  roughness"  of  the 
gas-liquid  ’nte-'face  At  a  critical  value  of  W*  the  inception  of  film 
instability  occurs,  and  a  further  increase  in  W*  causes  the  phenomenon 
of  coolant  entrainment  which  results  in  an  even  more  pronounced  reduction 
in  the  effectiveness  of  liquid-film  cooling  Because  of  these  points, 
and  since  the  heat  transfer  phenomena  characterizing  liquid-film 
cooling  and  transpi ration  cooling  would  appear  to  become  progressively 
more  similar  as  W+  is  decreased,  the  following  correlation  is  suggested 
for  experimenta1  heat-mass  transfer  data  that  are  obtained  for  liquid- 
fi lm  coo l mg  In  us 

>  ,(v/)  ]t  (3-?4> 

0  1  f  0 

limit  ffW+)  -  1 
W  •  o 


where 
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The  notation  and  subscripts  in  Eq.  3-24  denote  the  following: 

St  =  the  Stanton  number  for  heat  transfer 
f(W+'r  a  function  that  must  be  determined  experimentally. 

If  -  evaluated  fo^  the  case  of  liquid  film  cooling 

t  -  evaluated  for  the  case  of  transpi ration  cooling 

o  -  evaluated  for  the  zero-mass  transfer  case 

The  Stanton  number  for  heat  transfer'  is  defined  by  the  following 
expression  (29).  Thus 


St  - 


where 


<ou>e  cp 


the  recovery  temperature  for  the  main  stream, 
the  recovery  factor. 


(3-25) 


(3-26) 


As  shown  by  Spalding  (29),  the  recovery  factor  r  can  be  determined  in 
most  instances  w-'th  sufficient  accuracy  from  the  following  empirical 
equation  due  to  Squire.  Thus 


(3-27) 


The  subscript  notation  (  )*  introduced  in  Eq  3-27  means  that  the 

subscripted  quantity  is  to  be  evaluated  at  an  appropriate  reference  state. 
The  reference  state  suggested  herein  is  that  defined  by  the  reference 
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temperature  t*  due  to  Eckert  (30),  and  the  reference  concentration  C  * 
due  to  Knuth  (31).  Thus 


+  V  +  r 


(3-28) 


+  C  ) 
v,s' 


(3-29)1 


The  function  f(W+)  in  Eq.  3-24  must  be  determined  experimentally. 
Procedures  are  outlined  in  Section  3.3.4  for  evaluating  f(W+)  given 
either  (a)  bulk  heat  transfer  data  (i.e.,  m‘o  as  a  function  of  L  for 
various  main-stream  conditions  and  liquid-film  coolants),  or, 
preferably,  (b)  the  streamwise  distribution  of  (St)lf  as  determined 
from  detailed  measurements  of  the  structure  of  the  boundary-layer  region 
adjacent  to  the  surface  of  the  liquid  film. 

:o  facilitate  implementing  the  correlation  procedure  suggested  by 
Eq  3-24,  Section  3.3.1  presents  a  method  for  obtaining  an  approximate 
evaluation  of  the  temperature  at  the  gas-liquid  interface,  and  Sections 
3,3.2  and  3.3.3  suggest  methods  for  evaluating  (St/StQ)t  and  St0, 
respectively. 

It  should  be  noted  that  the  correlation  procedure  developed  herein 
differs  from  the  semi -empirical  correlations  reviewed  in  Section  3.1 
in  that  it  is  applicable  .o  (lows  in  the  presence  of  pressure  gradients, 
to  the  flow  of  a  compressible  fluid,  and  to  either  stable  or  unstable 

*  In  the  developments  which  follow,  it  is  assumed  that  the  concentration 

of  the  vapor  in  the  main  stream  is  zero.  Hence,  C  *  *  H  (Cw  .). 

V  )  v  »s 
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liquid  films  Moreover,  "entrance  length"  effects  and  streamwoe 
variations  in  q<.  and  { p v ) ^  are  considered.  The  developments  are 
limited,  however,  to  the  case  where  the  gas  flow  is  turbulent  and  the 
liquid-film  coolant  is  non-reactive. 

3.3.1  A  Method  for-  Predicting  the  Temperature  at  the 
Gas-Liquid  Interface 

The  temperature  at  the  gas-liquid  interface,  te  ,  can  be  related  to 
the  bulk  flow  parameters  by  specifying  (a)  an  interfacial  energy 
balance,  and  (b)  an  analogy  between  the  rates  of  heat  and  mass  transfer. 
In  the  derivations  that  follow,  the  generalized  energy  balance  given  by 
the  following  equation  is  employed  for  convenience: 

q<.  <t  Lv)s  (3-30) 

The  parameter  ®  is  defined  by  Eqs.  3-30  and  2-31.  Thus 


t 


:.H 

v 


-  qriS)/(°v)5 


The  analogy  between  heat  and  mass  transfer  employed  herein  is  that 
given  by 


St  =  St' 


(El)  ' 

'Sc'* 


2/3 


(3-31) 


where  St’  is  tne  Stanton  number  for  mass  transfer  defined  by 


St' 


Uv\ 


p-p 


V,S 


v,s 


R 


(3-32) 


1 


1 


Equation  ?-32  also  can  be  written  in  terms  of  the  concentration  of  the 
vapor  3$  f o 1 1 o vs  (29)  Thus 


(7uT 


Vj± 

v,s 


St'  £ 
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Experimental  verification  of  the  validity  of  Eq.  3-31  for  incompressible 
flows  is  presented  in  Ref  (32),  and  for  compressible  flows  in  Ref  (29). 
It  should  be  noted,  however,  that  at  and  close  to  the  plane  where  the 
liquid  coolant  is  injected,  Eq  3-31  is  strictly  applicable  only  if  the 
growths  of  both  the  thermal  bounaary  layer  and  the  partial-pressure 
(or  vapor-concentration)  boundary  layer  start  at  that  plane. 

Combining  Eqs  3-30,  3-31,  and  3-32  yields 


(PX) 

Sc  * 


-2/3 


(3-33) 


In  Eq  3-33  there  are  three  unknowns:  o,  p  „  and  t  Tabulated  data 

v,s  s 

presenting  the  saturation  pressure  of  the  liquid  coolant  as  a  function 
of  its  saturation  temperature  represents  a  second  independent  relation¬ 
ship  for  pw  s  and  t$  Thus,  by  specifying  $,  t$  (or  its  equivalent, 

p  )  can  be  related  to  the  bu]k  flow  parameters 
v  ,s 

Figure  3-2  presents,  for  an  air/water  system,  calculated  values 

of  toe  interface  temperature  t$  as  a  functior  of  the  static  pressure  p 

and  the  main-stream  temperature  t  *  The  flow  conditions  assumed  in 

the  calculations  are  indicated  in  the  figure  The  saturation  curve 

for  water  presented  in  the  figure  illustrates  the  approximate  error 

introduced  by  the  assumption  that  t,  is  equal  to  the  saturation  tempera- 

ture  for  the  liquid  coolant  und  r  the  prevailing  static  pressure,  p, 

that  assumption  having  been  enployed  by  Refs  (6)( 1G)( il)(22)(24) 

1  The  calculations  were  performed  by  employing  Eq.  3-33  in  conjunction 
with  the  tabulated  data  of  Refs  ( 3° )  and  (33)  A  sample  calculation 
is  presented  in  Appendix  B 
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TEMPERATURE  AT  INTERFACE  ,ts(R) 


TEMPERATURE  AT  INTERFACE  ,t6 (R) 


FIG.  3-2.  THE  TEMPERATURE  AT  THE  INTERFACE  AS  A 
FUNCTION  OF  THE  STATIC  PRESSURE  AND 
THE  MAIN-STREAM  TEMPERATURE 
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The  calculated  results  presented  in  Fig.  3.2  agree  within 
±  10F  with  the  data  reported  in  Ref.  (8)  and  that  obtained  in  the 
preliminary  investigations  of  liquid-film  cooling  conducted  at  the  Jet 
Propulsion  Center,  Purdue  University  (12);  the  flow  conditions  presented 
in  the  figure  apply  closely  to  those  experimental  investigations. 

The  values  of  p  and  tg  employed  in  the  subject  investigations  were 
as  follows: 

Ref.  (8):  t  =  1000-230CR;  p  =  1  atm 

Ref  (12):  t  =  1000k,  p  =  220  psia 

It  was  stated  in  Section  2.2.2  that  there  is  a  question  (especially 

for  the  case  where  a  streamwise  pressure  gradient  is  imposed  on  the  flow) 

regarding  the  fraction  of  the  energy  incident  to  the  surface  of  the  liquid 
film  that  is  utilized  for  increasing  the  sensible  enthalpy  of  the  liquid 
coolant  The  method  presented  herein  for  evaluating  t$  affords  a  means 
by  which  that  question  can  be  investigated  experimentally.  The  experi¬ 
mental  measurement  of  the  streamwise  distribution  of  t$  can  be  related 
to  a  streamwise  variation  in  the  parameter  by  means  of  Eq.  3-33. 
Moreover,  if  the  radiant  and  wall  heat  fluxes  are  reduced  in  importance, 
by  proper  design  of  the  experimental  apparatus  and  selection  of  the 
experimental  parameters,  then  the  definition  of  4  reduces  to 

<j>  -  aHv  ♦  (the  increase  in  the  sensible  enthalpy  for  the  liquid 
coolant) 

Thus,  for  that  case,  the  fraction  of  the  incident  heat  flux  utilized 
for  increasing  the  sensible  enthalpy  of  the  liquid  can  be  determined 
as  a  function  of  x 
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3.3.2  Correlation  of  ( St /StQ ) t  for  Transpiration  Cooling 

f 

Analytically  and/or  experimentally  it  has  been  established  satis¬ 
factorily  that  if  transpiration  cooling  data  are  presented  in  the  general 
form  1  -  — 


as  a  function  of 


'ov)<. 

(ou)e 


the  resulting  correlation  is  essentially  independent  of  the  Reynolds 
number  for  the  main  stream,  based  on  the  streamwise  coordinate  x,  the 
Mach  number  for  the  main  stream,  and  th6  ratio  of  the  main-stream 

i 

temperature  to  the  temperature  at  the  wall  (34} (35) (36).  Furthermore, 
the  correlation  is  independent  of  the  streamwise  pressure  gradient 
imposed  on  the  flow  (36),  Because  of  the  significance  of  those  facts, 
such  a  correlation  for  transpiration  cooling  data  is  employed  In  the 
present  developments.  j 

Figure  3-3  presents  the  correlation  for  transpiration  cooling  data 
due  to  Brunner  (37).  The  curve  presented  in  the  figure  is  given  by  the 
following  equation.  Thus 


(.— -) 

lSt  ' 


o  t 


<pv)s 


where 


0.6 


(3-34) 


I 


(3-35) 
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The  data  surveyed  by  Brunner  were  for  the  case  where  either  nitrogen, 
or  helium,  or  water  vapor  is  injected  into  an  uir  stream  As  a  con¬ 
sequence,  there  is  no  assurance  that  the  empirical  correlation  given 
by  Eq.  3-34  is  applicable  to  all  flowing  gases  and  injected  fluids  that 
might  characterize  a  transpiration-cooled  system  The  fluids  N^,  Hg, 
and  H^O  do,  however,  represent  a  significantly  wide  cross  section  of 
physical  property  values.  It  appears  reasonable,  therefore,  that 
Eq.  3-34  can  be  extended  with  some  degree  of  confidence  to  other  fluids 

333  A  Method  for  Evaluating  StQ 

IdeaMy,  in  an  investigation  of  liquid-film  cooling,  the  streamwise 
distribution  of  the  Stanton  number  for  heat  transfer  in  the  absence  of 
mass  transer,  denoted  by  Sto,  would  be  determined  experimentally.  There 
are,  however,  many  instances  where  experimental  determination  of  StQ  is 
impractical.  In  those  cases,  the  method  suggested  herein  can  be 
employed  for  ^yaluating  StQ  analytically 

The  method  suggested  has  been  developed  and/or  employed  by  numerous 
investigators  ( 36 ) ( 38 ) ( 39 ) ( 40 ) ( 4 1 ) ,  It  is  based  on  the  assumption  that 
the  empirical  heat- transfer  law  for  turbulent  flow  over  a  flat  plate, 
when  written  in  the  following  form,  is  also  applicable  to  turbulent  flows 
in  the  presence  of  streamwise  pressure  gradients.  Thus 

Sto  *  Sto  <V 

where 

p  u  i 

R  i  - — —  =  the  Reynolds  number  for  the  main  stream  based  on 

*  ^ 


the  characteristic  length 
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t  iou>  (tnu  -  t°; 

£  •  o  Ioule  (t~T-|gT  uy  ^3-36) 

=  tne  enerqy  tnicxness  . 

At=  the  thermal  boundary  layer  thickness  « 

2 

t°:  t  +  2~—  =  the  stagnation  temperature  . 

P 

The  heat-transfer  law  suggested  ir  Ref..  (36),  and  employed  in  the 
followino  deveiopments ,  is  that  oiven  bv 

StQ  =  0.0128  (Pr)‘3/4  R^“1/4  (3-37) 


tquation  3-37  and  the  developments  that  follow  can  be  employee  to 


analyze  a  compressible  qas  flow  by  introducing  the  concept  of  the  ref¬ 
erence  state,  and  making  the  following  transformations  for  the  Stanton 
number  StQ  and  the  Reynolds  number  Rf  (42).  Thus 


(3-38) 


(3-39) 


The  energy  integral  equation  is  given  by  (38) 


Ar  !T  dUQ 

Ct  3  t  r  l-i - §. 

blo  dx  f-iuD  dx 


d<VV  ,  1  - 

dx  p.  dx 


(3-40) 


Equation  3-37  substituted  into  Eq  3-40  results  in  an  explicit  equation 
for  StQ  that  can  be  solved  stepwise  by  numerical  integration  To  start 
the  numerical  calculations,  however,  the  value  of  xQ  and  the  value  of 
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StQ  at  x  =  xQ  must  be  either  known  (determined  experimentally)  or 

approximated  In  many  applications  of  liquid-film  cooling  that  are 

of  interest,  the  temperature  of  the  liquid  film  is  significantly  less 

than  the  temperature  of  the  solid  surface  upstream  of  the  plane  at  which 

film  coolant  is  injected-  In  those  cases,  it  is  suggested  that  StQ 

be  calculated  as  if  the  growth  of  the  thermal  boundary  layer  began 

at  x  --  xQ.  In  oth^r  words,  assume  that  the  boundary  condition 

St  =  ®  at  x  =  x„  applies, 
o  o 

In  the  special  case  of  an  incompressible  gas  flow  over  an  iso¬ 
thermal  wall  surface,  Eq.  3-40  can  be  simplified  to  the  form 


(3-41) 


Substituting  Eq  3-37  into  Eq  3-41,  integrating  the  resultant 
expression,  and  imposing  the  limit  StQ  =  ®  at  x  =  xQ  ;one  obtains 
the  following,  relatively  simple  expression  for  StQ.  Thus 


St 


o 


0.02923  Pr‘3/5  Rx"1/5 


v 


x 

u„  dx 
e 


(3-42) 


where 

°e  ue  x 

R  - -- - =  the  main-stream  Reynolds  number  based  on  the 

ue 

characteristic  length  x. 


Equation  3-42  reduces  to  the  following  expression  for  the  case 
where  the  main-stream  velocity  ug  is  constant  in  the  streamwise 


direction.  Thus 
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St0  =  0.02923  Pr'3/5  Rx'1/5  (1  -  —)  "1/5 


(3-43) 


x  -1/5 

The  term  (1  -  ~  )  in  Eq.  3-43  is  the  correction  for  "entrance 
length"  effects  due  to  the  dynamic  entrance  length  x  . 


3.3,4  Methodology 

The  principal  equations  required  for  the  correlation  procedure 
suggested  herein  are  summarized  below.  Thus,  by  definition 


U1 


(2-10) 


From  the  conservation  of  film  coolant  cne  obtains 


m'  =  m^  -  '  (Pv)s  dx 


o 


Combining  Eqs  3-25  arid  3-30,  one  obtains 


(ov)r 

Stlf  =  J 


(3-44) 


(3-45) 


whe  re 


J  5 


cp  (tr  ‘  ‘s’ 
Ke 


(3-46) 


Furthermore,  combining  Eqs  3-24,  3-34,  and  3-45,  yields 


-3 


iff) 

If 


f(w‘)  !  i  *  I  (If) 

I  J  Sto  If 


(3-47) 
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In  the  discussion  that  follows,  it  is  assumed  that  the  streamwise 
variations  for  the  quantities  listed  in  Table  2  are  either  determined 
experimentally,  or  can  be  approximated  from  knowledge  of  the  experiment¬ 
al  apparatus  and  parameters,  or  can  be  determined  analytically  from  the 
procedures  presented  in  Sections  3.3.1  and  3.3.3. 

Table  2 

The  Parameters  That  Are  Considered  Known  As  a  Function  of  x 
In  Outlining  the  Methodology 


Symbol  Definition 

c  the  specific  heat  for  the  main-stream  gas 

Pe 

c  the  specific  heat  for  the  vapor 

Pv 

I  non-dimensional  parameter  defined  by  Eq.  3-35 

J  non-dimensional  parameter  defined  by  Eq.  3-46 

L  the  liquid-film-cooled  length 

p  fhe  static  pressure 

St„  the  Stanton  number  for  heat  transfer  in  the  absence 

o 

of  mass  transfer 

t  the  main-stream  temperature 

ts  the  temperature  at  the  gas-liquid  interface 

u  the  velocity  of  the  main  stream 

e  J 

xQ  the  plane  at  which  liquid-film  coolant  is  introduced 

A  suggested  procedure  for  evaluating  the  function  f(W+)  is  outlined 
below  for  the  following  two  cases:  (a)  tne  case  where  the  streetwise 
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distribution  of  (St)^  is  determined  from  detailed  measurements  of 
the  structure  of  the  boundary- layer  region  above  the  liquid  film; 
and  (b)  the  case  where  only  bulk  heat  transfer  data  are  obtained 
(i.e.,  m'  as  a  function  of  L  for  various  main-stream  conditions 
and  liquid  film  coolants) 

Case  (a):  (St)^  Determined  Experimentally  The  procedure  for 

this  case  can  be  outlined  as  follows.  Thus 

(1)  Determine  the  streamwise  distribution  of  (p v ) ^  from  Eq  3-45 - 

(2)  Determine  the  streamwise  distribution  of  m'  and  W+  from  Eqs.  3-44 
and  2-10,  respectively 

(3)  Determine  the  streamwise  distribution  of  f(W+)  from  Eq.  3-47. 

From  steps  (1),  (2),  and  (3)  one  obtains,  for  each  set  of  input 

data,  a  plot  of  f(W+)  as  a  function  of  W+ .  Hence,  for  n  sets  of  input 
data  one  obtains  n  such  plots  If  the  basic  premise  upon  which  the 
correlation  suggested  herein  is  correct,  that  is,  if  f(W+)  is  an 
universal  function,  each  of  those  plots  would  be  identical  If, 
however,  the  curves  differed  significantly,  it  would  indicate  that 
the  function  f(W+)  is  not  sufficient  to  characterize  the  difference 
between  the  heat  transfer  phenomena  for  liquid-film  cooling  and  trans¬ 
piration  cooling  However,  even  if  that  was  found  to  be  the  case,  if 
f(W+)  differed  from  unity  in  any  instance,  the  inadequacy  of  employing 
solely  transpi ration  cooling  theory  in  the  analyses  of  liquid-film 
cooling  would  have  been  established 

Case  (b):  fiT  and  L  Determined  Experimentally  In  this  case  the 
functional  form  of  f(W+)  is  assumed;  i  e  ,  the  function  f(W+)  is  assumed 
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to  be  linear,  exponential,  etc  By  studying  several  sets  of  input  data, 
it  can  be  determined  whether  or  not  the  assumed  form  of  f(W+)  was 
satisfactory .  The  complete  procedure  is  outlined  below.  Thus 

(1)  Assume  the  functional  form  of  f(W+). 

(2)  Determine  m‘  and  W+  at  x  =  x 

o 

(3)  Calculate  (St)w  at  x  =  x  ,  employing  Eq.  3-47. 

i  <  o 

(4)  Calculate  (pv)<.  at  x  =  xq  by  means  of  Eq.  3-45. 

(5)  Determine  m'  and  W+  at  x  =  xQ  +  ax,  employing  Eqs.  3-44  and  2-10, 

respectively. 

(6)  Repeat  steps  (2)  through  (5),  replacing  x  =  xQ  by  x  =  xQ  +  ax, 

and  continue  the  reiteration  until  x  =  x,  =  x  +  L. 

1  o 

(7)  If  m' f  Oat  x  =  x^,  then  changes  must  be  made  in  step  (1).  To 

illustrate,  if  f(W+)  was  assumed  to  be  linear  in  step  (1),  so  that 

f(W+)  =  1+CW+,  then  a  different  value  must  be  assumed  for  the  constant 

C  and  the  computation  repea*°d  until  m'  =  0  when  x  =  x^. 

For  a  single  set  of  input  data,  regardless  of  the  assumed  func¬ 
tional  form  for  f(W+),  the  condition  m'  =  o  at  x  =  Xj  will  be 
realized  By  considering  several  sets  of  input  data,  however,  it 
can  be  deduced  whether  that  form  for  f(kT)  satisfies  the  condition 
of  universality  Note  that  if  one  experiences  difficulty  in  obtain¬ 
ing  a  form  for  f(W+)  tnat  simultaneously  satisfies  several  sets  of 
input  data,  it  is  nearly  impossible  to  determine  whether  that  difficulty 
is  due  to  a  lack  of  insight  into  the  physics  of  the  problem  (i.e.,  the 
inability  t  ’guess"  at  the  correct  f(W*)},  or  to  the  fact  that  f(W  ) 


does  not  physically  exist. 
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4.  AN  ANALYSIS  TOR  THE  GAS-VAPOR- COOLED  REGION 

4  1  Introduction 

It  is  difficult,  if  not  impossible,  to  make  a  complete  mathematical 
analysis  of  the  boundary-layer  and  heat-transfer  characteristics  for 
the  region  immediately  downstream  ^m  the  end  of  the  liquid  film. 
Moreover,  an  involved  analysis  generally  -is  not  warranted  because  of 
the  large  number  of  assumptions  that  must  be  made  concerning  tie  velocity 
and  temperature  profiles,  the  form  factor,  the  skin-friction  law,  the 
heat-transfer  law,  etc  As  a  consequence  of  this  fact,  the  analysis 
presented  herein  is  based  on  an  idealized  and  relatively  simple  flow 
model 

Both  Ennons  (6)  and  Guinn  (7)  developed  semi -empi ri cal  analyses 
for  the  gas- vapor-cooled  region;  those  analyses  were  based  on  the  flow 
model  suggested  by  Hatch  and  Papeii  (43)  for  gas-film  cooling  Tne 
analysis  presented  here  differs  from  those  due  to  Emmons  and  to  Guinn 
in  the  following  respects:  (a)  bounda  ry- 1  ayer  flow  rather  tnan  fully- 
developeu-pi pe  flow  's  considered;  and  (b)  the  analysis  is  generalized 
so  that  it  can  be  applied  to  a  compressible  gas  flow  in  tne  presence 
of  a  streanwise  pressure  qradient  The  analysis  is  dependent,  however, 
on  tne  analysis  for  the  1 iquid-fi Im-cooled  region  presented  in 


Section  3.3. 
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As  discussed  in  Section  2  12,  a  complete  analysis  for  the  gas-vapor- 
cooled  region  requires  the  consideration  of  two  cases:  (a)  the  case 
where  the  wall  is  adiabatic;  and  (b)  the  case  where  the  wall  is  externally 
cooled.  Those  two  cases  are  considered  in  Sections  4.2  and  4.3, 
respecti  vel  y . 


+  ]iie  Case  Where  the  Wall  js_Adiabatic 
4.2.1  The  Flow  Model 

The  flow  model  employed  for  the  analysis  of  the  gas- vapor-cooled 
region  is  defined  l>v  the  fol  lowina  conditions 

( 1)  Tne  wal 1  is  adiabatic. 

(2)  Tne  boundary  layer  flow  is  fully  turbulent 

(3)  The  thermal,  tie  diffusional,  and  the  dynamic  boundary  layer 
tni tkrusses  are  equal,  the  latter  is  denotea  by  y  =  A. 

{•;)  Tne  stagflation  enthalpy  for  tne  main  stream,  denoted  by  h^,  is  a 
constant  in  tne  streamwise  direction;  the  stagnation  enthalpy  for  the 
qas  stream  is  defined,  in  general,  by 

y  *  f  <«-» 

(5)  The  velocity  profile  and  the  stagnation-enthalpy  profile  for  the 
boundary  layer  are  of  the  functional  form’s 


g  .  * 


ana 


a*. 


■  f,  (i ) 


h  -  n 


aw 


(4-2) 


(4-3) 
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where  h  ’s  ti  “  or  thalpy  for  tie  fluid  adjacent  to  the  adiabatic 
wall.  Tc  illustrate,  f^  could  be  assumed  of  the  power-law  form: 

1 

f!  *  If)"  (4-4) 

In  addition,  could  be  assumed  of  the  form 

U  =  A  ♦  B  (*)  ♦  C  (^)2  ♦  D  (*)  (4-5) 

&•  ci  Q 

The  constants  in  £q .  4-5  can  be  evaluated  from  the  following  boundary 
conditions  that  apply  to  f~.  Thus 

f2  (0)  *  0 

f2  (1)  1 

hf, 

if{0)  =  0 

;f, 

~(1)  ’  0 

I’  ■  ultant  expression  for  f  is 
i  3 

M\)  -  2{*)  (4-6) 

Tor  tfie  present  analysis,  however,  the  actual  forms  of  f,  and  f„ 

1  c 

nee^  not  be  assumed,  only  tne  assumption  tnat  f.  and  f?  exist  for  the 

1  K 

gas- vapor-cooled  re/>cn  is  required. 
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4  2.2  The  Energy  Balance 

Figure  4-1  presents,  schematically,  a  control  volume  for  the  gas- 
vapor-cooled  region;  the  control  surface  for  that  control  volume  is 
defined  by  x  -  ,  x  -  x,  y  =  ,  an j  y  =  o.  The  foliowinq  energy 

fluxes,  written  in  accordance  with  the  flow  model  defined  in  the  fore¬ 

going  section,  are  presented  in  the  figure: 

£ 

1 

(a)  eu  h°  dy  =  the  rate  at  whic:  thermal  and  kinetic  energy  is 

o 

transferred  into  the  control  volume  per  unit 
spanwise  length  at  the  station  x  =  x^. 

(b)  cu  h°  dy  -  the  rate  at  which  thermal  and  kinetic  energy  is 

0  transferred  out  of  the  control  volume  per  unit 

spanwise  length  at  the  station  x  =  x 

(c)  (m"  -  m'j)h°  -  the  ^te  at  which  thermal  and  kinetic  energy  is 

transferred  from  the  main  stream  into  the  control 
volume  per  unit  spanwise  length 
The  parameter  m“  in  term  (c)  is  defined  as  fellows: 

m“  cu  dy  (4-7) 


-  trie  rate  at  which  mass  is  transferred  throuqn  the  bouncer* 
layer  pe-  unit  spanwise  length  at  fie  station  x  - 
The  energy  balance  for  the  control  volume  in  Fig  4-1  can  l«- 
written  as  follows  Thus 


•  ^  0 

cu  n  Gy  ♦  (r  -r  ^  •- 


c  j  h  dy 


(4-6) 


STREAM 
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FIG.  4-1.  THE  CONTROL  VOLUME  FOR  THE  GAS-VAPOR-COOLED 
REGION 
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Employing  Eq,  4-7,  Eq.  4-8  can  be  written  in  the  form 


'  .0  . 
i  pu  h  dy 


/  pu  dy 


+  h°  =  m" 


pu  h°  dy 


/  pu  dy 


(4-9) 


In  order  to  realize  the  desired  simplicity  of  the  present  analysis,  the 

following  approximation  is  introduced  concerning  the  integral  terms  in 

Eq,  4-9,  Thus* 

A  A 

/  pu  h  dy  /  u  h°  dy 

0 _  a  0 _ 

A  '  A 

/  pu  dy  /  u  dy 


Thus,  Eq.  4-9  can  be  written  in  the  following  form  by  introducing  that 
approximation  and  Assumption  (5)  of  the  flow  model: 


•  V(he  -  haw.l>Wl> 

^ _ a  \  lO  _ 


+  (m"-m1‘)  h”  = 


/  dy 


/  fl  «he  -  haw>Vhaw>  « 

O _ 

A 

/  fx  dy 


The  latter  expression  can  be  rewritten  in  the  form 


Note  that  this  approximation  is  not  equivalent  to  and  is  less 
restrictive  than  the  assumption  of  constant  density., 
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m"  {h°  -  h  ,)  -  m"  (h°  -  h  ■  .)  - 
1  v  e  aw,l'  1  '  e  aw,l' 


/  fj  f2  <*(£> 


l  'l 


m"  (h°  -  h  ) 
e  aw 


m"  (h°  -  h  ) 

'  e  aw'  i 


i 

0  1  c 


'fl  d<^) 


(4-10) 


Finally,  the  terms  in  Eq  4-10  can  be  combined  to  yield  the  following 
energy  balance  for  the  control  volume  presented  in  Fig,  4-1.  Thus 


h 


h 


e 

o 

e 


m 


m' 


(4-11) 


The  relatively  simple  form  of  Eq  4-11  is  due  principally  to  Assumption 
(5)  of  the  fnow  model , 

Equation  4-11  enables  the  determination  of  the  enthalpy  for  the 
fluid  adjacent  to  the  adiabatic  wall  as  a  function  of  x  once  h°  and  h 

“  dr 

and  nf'  as  functions  of  x  are  known  The  parameter  h°  is  normally 
known  from  the  statement  of  the  problem,  and  the  parameter  h  .  can  be 

oW  j  1 

determined  from  the  developments  presented  in  Section  3  3  1  The 
evaluation  of  m"  as  a  function  of  x  is  considered  in  Section  4  2  3 
To  evaluate  the  temperature  distribution  for  the  adiabatic  wall 
from  the  distribution  of  h  ,  the  specific  heat  for  the  fluid  adjacent 

dw 

to  the  wall  must  be  known.  Excluding  the  region  close  to  the  station 
x  =  Xj,  it  is  sufficiently  accurate  to  take  that  specific  heat  to  be 


bb 

that  for  the  main-stream  gas  because  of  the  low  concentration  of  the 
vapor  in  comparison  to  that  of  the  main-stream  gas.  Moreover,  since 
the  temperature  of  the  adiabatic  wall  at  x  =  x^  can  be  independently 
determined  from  the  developments  of  Section  3.3.1,  it  follows  that  the 
adiabatic-wall  temperature  distribution  for  the  entire  gas- vapor-cooled 
region  can  be  determined  with  sufficient  accuracy  once  li  as  a  function 
of  x  is  known. 


4.2.3  The  Evaluation  of  m" 

Equation  4-5,  which  defines  the  parameter  m" ,  can  be  rewritten  in 
the  form 

m"  *  <pu)e  <1  •  r>  0  <4-12> 

where 

a 

a*  z  j  (1  -  -^y  )dy  =  the  displacement  thickness.  (4-13) 

o  e 

A 

0  =  /  T^-y (1  -  j-j-Jdy  =  the  momentum  thickness.  (4-14) 

o  p  'e  e 

Equation  4-12  shows  that  m"  can  be  determined  as  a  function  of  x  if 
A/e,  A*/0.  and  e  are  known  as  functions  of  x.  The  evaluation  of  those 
quantities  is  considered  in  the  following  paragraphs. 

The  Evaluation  of  A/e  and  A*/e.  Persh  and  Lee  (44)  have  tabulated 
the  values  of  A/e  and  a*/6  as  functions  of  the  main-stream  Mach  number 
and  the  wal 1-to-main-stream  temperature  ratio.  It  is  suggested,  there¬ 
fore,  that  these  tabulated  data  be  employed  for  evaluating  A/e  and 
A*/e  when  the  variation  in  the  density  for  the  gas  stream  across  the 


i 
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boundary  layer  is  significant.  If,  however,  that  density  variation  is 
not  significant,  it  is  suggested  that  the  subject  parameters  be  calculated 
by  applying  the  power-law  velocity  profile  (see  Eq.  4-4).  The  resultant 
expressions  for  A/e  and  a*/0  are  (4,  p.  536) 


_  (2+n) 
e  n 


(4-15) 

(4-16) 


The  Evaluation  of  e.  To  determine  the  value  of  the  momentum 
thickness  e  as  a  function  of  x,  one  must  solve  the  momentum  integral 
equation.  That  equation  can  be  written  in  the  form  (45) 


c_f  _  de 
2  ~  dx 


♦  e  [(*♦£)  L  Zl  ♦  JL 
u  e  '  u  dx  o 
e  e 


(4-17) 


Or  in  the  form  (46) 


f  w 

__  - -  =  the  skin-friction  coefficient.  (4-19) 

De  ue 

M  =  the  Mach  number  for  the  main-stream  flow, 
c 

\  i  -  the  specific  heat  ratio. 
cv 
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To  solve*  either  Eq.  4-17  or  Eq,  4-18  for  o,  one  muse  specify 
the  form  of  the  skin-friction  law;  the  latter  is  denoted  functionally 
by 


where 

p  u  e 

d  5  - -  -  the  Reynolds  number  for  the  mam  stream  based 

6  ue 

on  the  characteristic  length  e.  (4-20) 


Furthermore,  the  boundary  condition  on  8  at  the  station  x  *  x^, 
denoted  by  e^,  also  must  be  specified  before  e  as  a  function  of  x 
can  be  determined  from  either  Eq  4-17  or  Fq.  4-18 

Several  empirical  skin-friction  1  -ws  are  presented  in  Chapter  22 
of  Ref,  (4)  for  the  case  where  the  gas  flow  is  turbulent,  incompressible, 
and  over  a  flat  plate.  In  the  analysis  of  gas  flows  with  pressure 
gradients  (particularly  favorable  pressure  gradients),  it  is  common 
to  assume  that  the  skin-friction  law  is  that  established  empirically 
for  the  flat-plate  case  (45) (46) (4).  Furthermore,  the  case  where  the 
gas  flow  is  compressible  can  be  analyzed  by  introducing  the  following 
transformations  into  the  form  of  the  skin-friction  law  for  incompressible 
flow  (42): 


and 


A  numerical  solution  is  generally  required  for  either  Eq.  4-17  or 
4-18 
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R 


e 


0* 


(4-22) 


The  reference  state  denoted  by  (  )*  in  Eqs.  4-21  and  4-22  was  discussed 
in  Section  3.3.  It  is  suggested  that,  in  evaluating  the  physical 
properties,  the  composition  of  the  gas  be  assumed  uniform  throughout 
the  boundary  layer  and  equal  to  that  for  the  main  stream;  i.e  ,  assume 
that  the  reference  temperature  given  by  Eq  3-28  is  sufficient  to  define 
the  reference  state 

A  sufficiently  accurate  approximation  to  can  be  obtained  in 
the  following  manner  First,  assume  that  the  Colburn  analogy,  given 
by  the  following  equation,  applies  at  x  =  Thus 
cf  -2/3 

(St),f  *  f  (Pr),  M-23) 

Cf 

Then,  employing  the  skin-friction  ^w  which  relates  — ■  and  e,  one  can 
approximate  the  value  of  6j  once  (St)lf  at  x  =  x^  has  been  determined 
from  the  analysis  presented  in  Section  3  3  for  the  1  iquid-fi  lm-coolec4 
region 


424  A  Special  Case 

To  illustrate  the  developments  presented  in  the  foregoing,  consider 
the  specia'  case  defined  by  the  fo1 lowing  conditions 

1  The  gas  flow  is  incompressible,  constant-density,  and  over  an 
adiabatic  flat  plate 

2  The  skin-friction  law  is  that  due  to  Prandtl  (4,  Eq  22  7).  Thus 

cf 

2  =  0.012b  Rq 


(4-24) 
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The  velocity  profile  for  the  boundary  layer  is  given  by 
1/7 

-u-  =  (*, 
ue  V 


(4-25) 


4  The  specific  heat  for  the  vapor  and  for  the  main-stream  gas 
are  equal  and  are  not  functions  of  temperature 

For  that  special  case,  Eq  4-11  reduces  to  the  form 


te~taw 

te'taw,l 


m. 


m 


(4-26) 


Equation  4-12  reduces  to  (refer  to  Eqs  4-15  and  4-16) 


m"  =  9  pue  e 


(4-27) 


Equation  4-17  simplifies  to  the  form 
c,  dR, 


’f 

~1 


de 

dx 


e 


dR 


(4-28) 


where 


ou  x 

Rx  i  — e  =  the  Reynolds  number  for  the  main  stream  based 

on  the  characteristic  length  x  (4-29) 

Upon  combining  Eqs  4-26  and  4-27,  one  obtains 


tetaw 
*e  ^aw.l 


(4-30) 


In  addition,  substitution  of  Eq.  4-24  into  Eq.  4-28  yields 


0  0128  dR  Ra4  d  Rc 

Xu  u 


(4-31) 


Integrating  Eq.  4-31  and  substituting  the  limit  R  ‘R.  at 

b  b  | 

x  -  Xj ,  one  obtains 
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Re5/4  -  0.016  Rx  (1  -  —)  +  R^/4 


(4-32) 


It  follows  from  Eq.  4-30  and  Eq.  4-32  that  for  the  special  case  in 
question,  the  adiabatic-wail  temperature  as  a  function  of  x  is  given 
in  non-dimensional  form  by  the  following  equation.  Thus 


te'taw  „  el 

te’taw,l  [0.016  Rx(l-£l)  ♦  R0i5/4]  4/5 


(4-33) 


Finally,  to  complete  the  solution,  RQ  can  be  related  to  (St)-p  p 
the  Stanton  number  for  heat  transfer  at  x  =  Xp  by  combining  Eqs.  4-21 


and  4-24.  The  resultant  expression  is 


-<st),fa,  W  -\ 


0  0128 


(4-34) 


4 . 3  The  Case  .Where  the  Wall  is  Cooled  Externally 
For  the  subject  case,  the  basic  analytical  problem  is  that  of 
relating  the  local  values  for  the  wall  temperature  and  the  wall-heat 
flux  to  the  main-stream  conditions  It  can  be  inferred  from  several 
mass-transfer  cooling  (and  heating)  studies  reported  in  the  literature 
(Refs.  47-51)  that  heat-transfer  data  that  are  obtained  for  the  gas- 
vapor-cooled  region  should  correlate  satisfactorily  in  the  following 
manner  Thus 

t  -  t 

St  =  St  — — - — 

gv  o  to  _  t 

e  w 


(4-35) 
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where 


St 


'w 


9«  i“9)e  (taw-tw) 


(4-36) 


-  the  Stanton  number  for  heat  transfer  for  the  gas-vapor- 
cooled  region 

StQ  =  the  Stanton  number  for  heat  transfer  In  the  absence  of 
gas- vapor  cooling 

-  the  temperature  the  wall  would  obtain  if  it  was 
adiabatic 

t  1  the  temperature  the  wall  obtains  due  to  the  fact  that  it 
is  cooled  external ly 

qw  -  the  heat  Hux  from  the  hot  gas  stream  to  the  external 
coolant 

Equation  4-35,  together  with  the  developments  presented  in 
Sections  4-?  and  3  3  3,  which  enable  the  determination  of  t...  and 
St^  as  functions  of  x,  respectively,  thus  represents  the  solution  for 
the  subject  case 


1  It  is  important  to  note  that  in  the  evaluation  of  St  in  Eq„  4-35, 
the  boundary  condition  StQ  -  ®  at  x  *  Xj  applies  (47551). 
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5  CONCLUSIONS  AND  RECOMMENDATIONS 

Having  studied  the  physical  and  theoretical  nature  of  liquid-film 
cooing,  together  with  the  pertinent  literature,  the  following  general 
conclusions  are  drawn: 

1  The  feasibility  of  liquid-film  cooling  a  rocket  thrust  chamber  has 
been  established  experimentally  (5)(6V(8)( 10) ( 11) (52 ' (53) (S4) . 

However,  the  fundamental  phenomena  characteri  zi  .fg  tnat  p.  icess  are  not 
well  understood.  In  particular,  the  mterfac'al  ,  lenonier  (film 
instability,  droplet  entrainment,  and  the  "effective1'  film  roughness) 

j  not  well  defined  especially  for  the  case  where  the  pressure  and/or 
main-stream  temperature  ij  substantially  gr  >ater  than  ambient  It  is 
generally  agreed,  however,  that  those  phenomena  significantly  influence 
the  effectiveness  of  liquid-film  coolin 

2  A  need  exists  for  accurate  heat-mass  transfer  data  for  the  liquid- 
nlm-cooled  region,  especially  fo-  the  case  characterized  by  one  or  more 
of  the  to’ ’owing  conditions:  (a)  the  ga^  flow  is  compressible;  (b)  a 
pressure  gradient  is  imposed  on  the  flow;  and  (c)  a  film  coolant  otn»r 
tnan  water  is  employes 

3  The  erases  presentee  to  date  for  liquid-film  cooling  have  either 
lacked  agreement  w’th  the  experimental  data  or  have  resulted  m  correla¬ 
tions  that  are  net  justified  physically  obviously,  further  experimental 
evaluation  of  tne  basic  pnenomena  would  enable  the  postulation  of  more 


93 


realistic  analytical  models  and  the  develop. .ent  of  mere  complete 
analyses 

4  Additional  experimental  investigation  of  the  gas-vapor-cooled 
region  is  needed  Boundary  layer  measurements  would  aid  in  defining 
the  mixing  characteristics  for  that  region,  and  would  enable  the 
postulation  of  a  more  realistic  f<ow  model  than  those  postulated 
previously  (6)(7),  and  that  presented  in  the  subject  report 

In  view  of  the  aforementioned  conclusions,  and  in  accordance 
with  the  foregoing  developments  presented  in  the  subject  report, 
some  recommendations  for  further  research  are  presented  below: 

1.  Conduct  cold-flow  studies  <  i.e  ,  zero  evaporative  mass  trans^r) 
of  film  instability  and  entrainment  to  determine  the  influence,  if 
any,  of  the  dimensionless  liquid  flow  rate  W*,  the  dimensionless  film 
thickness  the  Reynolds  number  for  the  gas  stream,  the  pressure, 
etc  on  those  phenomena  For  these  studies,  and  for  those 
recommended  below,  several  fluids  in  addition  t"’  water  should  be 
studied,  and  pressures  significantly  greater  than  ambient  should 
be  employed 

l  Conduct  additional  co'd-flow  studies  utilizing  a  pitet  probe  to 
determine  the  velocity  profile  across  the  boundary  layer  cS  a  function 
of  x  Use  these  data  to  study  the  “effective1'  roughness  of  the  gas- 
liquid  interface 

3  Conduct  heat-mass  transfer  studies  to  obtain  data  that  can  be 
employee  to  evaluate  the  correlation  procedure  suggested  in  Section  3  3 
Initial  studies  should  be  conducted  to  obtain  bulk  heat-mass  transfer 
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data  for  a  significantly  wide  range  of  experimental  parameters. 
Subsequently ,  more  sophisticated  studies  should  be  conducted  to 
obtain  boundary  layer  data  that  would  aid  in  the  understanding 
of  the  phenomena  occurring  at  the  gas-liquid  interface 
4  Experimentally  investigate  the  region  downstream  of  the  liquid 
film  for  the  case  where  the  wall  is  adiabatic  Initial  studies 
should  be  conducted  wherein  only  the  wall  temperature  distribution 
and  the  bulk  flow  parameters  are  determined  Thor.e  data  can  be 
employed  to  evaluate  the  analysis  presented  in  Section  4,2 
Subsequently,  experimental  studies  should  be  conducted  to  obtain 
boundary  layer  measurements  that,  would  aid  in  the  postulation  of  a 
realistic  flow  model  for  the  subject,  wa1!  region 
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APPENDIX  A 

NOTATION 

specific  concentration  of  "i"  species 

skin  friction  coefficient  {Eq.  4-19) 

specific  heat  at  constant  pressure 

specific  heat  at  constant  volume 

tube  diameter 

diffusion  coefficient 

denotes  a  function 

specific  enthalpy 

heat  transfer  coefficient 

stagnation  enthalpy  (Eq  4-1) 

^tent  heat  of  vaporization 
e  Equation  3-35 
jee  -quation  3-46 
thermal  conductivity 
liquid-film-cooled  length 

rate  of  liquid  coolant  fiow/unrt  spanwise  length  (Eq  2-7) 
rate  of  gas  flow  through  the  boundary  layer/unit  spanwise 
length  (Eq  4-7) 
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Mach  number 

parameter  in  the  power- law  velocity  profile  (see  Eq.  4-4) 
static  pressure 

partial  pressure  for  "i"  species 

u  c 

Prandtl  number  =  -y-P- 

yt  CD 

turbulent  Prandtl  number  =  — r--*- 

t 

energy  transfer  due  to  conduction 
energy  transfer  due  to  radiation 
recovery  factor  (Eq.  3-26) 
perfect  gas  constant 


Reynolds  number  based  on  d 
Reynolds  number  based  on  x 


pe  ue  x 


Reynolds  number  based  on  £  =  - 


pe  ue  5 


Reynolds  number  based  or  9  = 


p„  LL  9 

e  e 


Schmidt  number 


turbulent  Schmidt  number  = 


Stanton  number  for  heat  transfer  (Eq.  3-2S) 

St  evaluated  in  the  absence  of  mass  transfer 
Stanton  number  for  mass  transfer  (Eq.  3-32) 
static  temperature 

recovery  temperature  for  the  main  stream  (Eq.  3-26) 


stagnation  temperature 
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\  .  saturation  temperature  for  the  liquid  coolant  which  cor- 

b  tl  L 

responds  to  the  static  pressure  p 
u  streamwise  component  of  velocity 

u+  dimensionless  velocity  (Eq.  2-5) 

v  normal  component  of  velocity 

W+  dimensionless  liquid  flow  rare  (Eq.  2-10) 

W+  see  Equation  2-37 

distance  downstream  from  the  effective  leading  edge  of  the 
dynamic  boundary  layer 

y  distance  normal  to  the  wall  or  the  liquid  film  surface 

y+  dimensionless  normal  coordinate  (tq.  2-b) 

Greek 

A  dynamic  Doundary  layer  thickness 

■\  diffusional  boundary  layer  thickness 

thermal  boundary  layer  thickness 

A*  disp iace-ient  thickness  (Eq.  4-13) 

heat  transfer  parameter  (Eq.  3-1) 

speci fi c  heat  ratio 

5  viscous  sublayer  thickness 

5 ^  liquid  fi 1m  thickress 

i*  dimensionless  viscous  sublayer  thickness 

* 

dimensionless  liquid  film  thickness  (Eq.  2-11) 

»■  stabi’ity  effectiveness  (Eq  2-38) 

J 

•  momentum  thickness  (Eq.  4-14) 


dynamic  viscosity 
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v  kinematic  viscosity 

S  energy  thickness  (Eq.  3-36) 

o  density 

t  shear  stress 

Ts  o  wa^  shear  stress  for  zero  mass  transfer 
$  energy  balance  parameter  (Eq.  3-30) 

Subscripts 

o  evaluated  at  the  plane  of  liquid  coolant  injection 

1  evaluated  at  the  terminus  of  the  liauid-film-coolec  length 

l  evaluated  it  the  terminus  of  the  gas-vapor-cooled  length 

aw  adi abati c  wal 1 

d  diffusional  component 

e  main-stream  state 

g  main-stream  gas 

yv  gas-vaDor  cooling 

1  liquid 

L  evaluatea  at  the  gas-liquid  interface  but  in  the  liquid  phase 

If  1  iquid-fi  1m  cool  ing 

s  evaluated  at  the  gas-liquid  interface 

evaluated  at  the  gas-liquid  interface  but  in  the  gas  phase 
t  transpi rati  on  coolinq 

t  turbulent  contribution 

v  vapor 

w  wall 


jvaulated  at  y  =  * 


105 


evaluated  at  the  reference  state  defined  by  Eqs.  3-28  and 
3-29 


Superscripts 

o  stagnation  conditions 
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APPENDIX  B 

A  SAMPLE  CALCULATION  FOR  THE  INTERFACE  TEMFERATURE 


To  illustrate  the  manner  in  which  the  calculated  results  t resented 
in  Fig.  3-2  were  obtained,  a  sample  calculation  tor  the  interface 
temperature  t.  is  presented  below. 

For  the  flow  conditions  presented  in  Fig.  3-2,  Eq,  3-33  reduces 
to  the  form 


AH 


V 


R.,  c-p 
=  _± 

P.  1  - 


) 


v,s 


83.35 
‘  53. 3!? 


p-p 


v,s 


V,s 


-  1  552  ( 


P-P 


v  ,s 


(B-l) 


Solving  Eq.  B-l  for  the  partial  press  ore  of  the  vapor  at  the  interface, 
one  obtains 


P 


v,s 


1 

1.562 


(B-2) 


For  purposes  of  illustration,  assume  that 

t  --  2000  R  (c 
e  d  ■ 

Kai  r 

p  =  500  psia 


=  0  2772  B/lb) 
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For  tia.  case,  Fq.  B-2  reduces  to 


500 

’v,s '  "T?r242TT3y  “ 

(2coo-TWT  +  1 

s 


psia 


(B-5) 


To  initiate  the  tr  al-and-error  calculations,  assume  that  t$  =  400F  =  860R 
(the  accuracy  of  the  first  approximation  for  t  is  not  critical).  For 
t  960°R,  the  data  in  Ref.  33  gives  AH  =  826.0  B/lb  and  p  =  247.3 

S  V  V )  s 

psia.  Substituting  those  data  into  Eq.  B-3,  one  obtains 


247. 3  ps «a  = 


500 


1  624  +  1 


psi  a 


/  190.5  psia 

For  the  second  calculation,  assume  that  p  =  190.5  psia.  Corresponding 

V  ,5 

to  p  *-  190.5  psia,  the  data  in  Ref.  33  gives  AH  =  846.7  B/lb  and 

V  yj  V 

t  :  837. 6R.  Substituting  those  data  into  Eq.  B-3,  one  obtains 

lnn  r  -  ?  500 

190  5  psia  =  \ —  Psia 

/  189  9  psia 


Continuing  that  procedure,  assume  for  the  next  calculation  that  p  = 

V  ,s 

189.9  psia.  Corresponding  to  p  =  189.9  psia,  the  data  in  Ref.  33 

V  *5 

gives  AH  =  846.9  B/lb  and  t  =  377. 4F  =  837. 4R.  Substituting  those 
data  into  Eq  B-3,  one  obtains 
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=  189.9  psia 

I 

Thus,  for  t  =  2000  R  and  p  =  500  psia,  and  for  the  flow  conditions  I 

S  j 

presented  in  Fig.  3-2,  the  calculated  value  of  the  temperature  at  the 

gas-liquid  interface  is  t$  =  377. 4F  =  837. 4R.  j 

I 

i 

I 

i 

I 

I 

! 

I 
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the  wall  is  cooled  externally. 
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